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I.  INTRODUCTION 

The  main  problem  of  sound  ranging  is  to  determine  the  location  of  a  tran¬ 
sient  sound  source  by  examining  the  signals  received  at  an  array  of  microphones 
The  problem  can  be  logically  broken  down  into  two  main  parts.  First  the  deter¬ 
mination  of  the  relative  arrival  times  of  the  sound  at  the  microphones  and 
second,  the  use  of  this  information  to  estimate  the  position.  in  this  paper 
we  concentrate  on  the  determination  of  relative  arrival  times. 

In  the  past,  the  arrival  times  at  each  microphone  were  determined  by  visu¬ 
ally  selecting  a  point  of  charge,  the  so-called  "break  point",  from  strip  chart 
recorded  signals  [1].  Some  methods  used  to  determine  these  break  points  have 
been:  first  inflection,  first  maximum  after  inflection  and  first  cross  over 
after  inflection.  A  comparative  study  of  these  methods  is  given  in  Dean  [2], 

In  contrast  to  the  break  point  methods,  the  author  [3]  proposed  a  procedure 
to  obtain  more  accurate  timing  information  from  the  received  microphone  sig¬ 
nals  by  using  the  correlation  between  signals.  The  method  when  implemented 
on  a  digital  computer  would  not  only  speed  up  the  determination  of  the  rela¬ 
tive  arrival  times  but  provide  more  accurate  results  with  less  chance  for  read¬ 
ing  errors  than  the  visual  techniques. 

The  purpose  of  this  paper  is  to  document  a  Fortran  realization  of  the  cor¬ 
relation  technique  described  in  [3]  and  present  results  from  using  field  data 
from  the  PASS  [4]  experiment  as  input.  These  limited  results  indicate  that, 
under  a  controlled  experimental  environment,  the  overall  method  provides  reason 
ably  good  time  differences  estimates  which,  when  used  with  the  proper  meteoro¬ 
logical  data,  give  good  position  estimates. 
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II.  GENERAL  DESCRIPTION  OF  TIME  DIFFERENCES  ESTIMATOR  PROGRAM 
The  overall  structure  of  the  computer  program  for  determining  the  rela¬ 
tive  time  differences  consists  of  the  following  basic  sections:  (A)  Initiali¬ 
zation  and  acceptance  of  input  data,  (B)  a  determination  of  the  time  windows 
that  would  enclose  the  signal  received  from  the  source  and  a  calculation  of 
the  energies  within  the  windows,  (C)  a  prefiltering  operation  that  would  over¬ 
come  wind  and  extraneous  noise,  (D)  a  pair  by  pair  correlation  to  determine  a 
rough  time  difference  and  a  measure  of  correlation,  (E)  the  determination  of 
an  overall  least  square  time  fit,  and  (F)  an  outputting  of  re:  'ts.  A  general 
flow  chart  indicating  processing  structure  is  given  in  Figure  1,  the  program 
listing  is  given  in  Appendix  A,  and  each  main  block  is  described  in  detail  in 
the  following  sections.  Each  of  the  subheadings  under  the  main  blocks  are  set 
apart  bv  comment  cards  in  the  program  for  ease  in  presentation. 

A .  Initialization  and  Acceptance  of  Data 

1.  Initializing.  In  the  initialization  portion  of  the  program  there  are 
a  number  of  parameters  that  must  be  specified  and  dimensioned.  The  main  param¬ 
eters  are: 

NMIC  =  Number  of  microphones 

NL  =  Number  of  data  points  for  each  microphone 
NSW  =  Number  of  points  in  the  sliding  window 
NPFW  =  number  of  points  in  the  Fourier  window 

The  number  of  points  in  the  sliding  window  is  chosen  to  be  approximately 
the  number  of  points  in  the  expected  signal.  By  examining  the  C4  data  from 
the  PASS  experiment,  an  example  of  which  appears  in  Appendix  C,  it  was  seen 
that  signals  were  in  the  300-375  msec  length  and  at  a  sample  rate  of  1  K 
samples/second  this  means  approximately  350  samples.  The  number  of  points  in 
the  Fourier  window  was  selected  to  be  1.5  times  the  number  in  the  sliding 


MICROPHONE  SIONAI.S 


INFORMATION  TIMES 

CORRECTED  TIME 
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Figure  1.  Processing  structure  of  the  Time  Differences  Estimator  Program 


window  to  make  sure  that  most  of  the  signal  was  within  the  window.  For  the 
C4  this  meant  approximately  500,  and  si-ace  the  Fourier  transform  can  be  ob¬ 
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tained  faster  for  powers  of  2  by  the  FFT,  NPFW  was  selected  to  be  512  and 
NWS  =  342.  The  data  files  provided  by  personnel  from  the  Physical  Sciences 
Laboratory  (PSL) ,  Las  Cruces,  New  Mexico,  were  constructed  to  be  1920  samples 
per  signal  in  length;  therefore,  ML  is  set  equal  to  1920.  The  following 
dimension  and  complex  cards  would  then  be  required: 

DIMENSION  IDAY(NMIC),  IHOUR(NMIC) ,  IMIN(NMIC) ,  SEC(NMIC) ,  MIC(NMIC) ,  TIME(NMIC) 

DIMENSION  TAU(NMIC,NMIC) ,  IT(NMIC) ,  NR(NMIC) ,  E(NMIC) ,  RTAU(NMIC,NMIC) , 

CTAU  (NMIC ,  NMIC)  ,  C.AM(NMIC  ,NMIC) 

DIMENSION  SMIC (NMIC, NL.LRST (NMIC,  G (NMIC, NMIC ) ,  AT (NMIC) ,  A (NMIC , NMIC) , 
JC(NMIC-l),  V(2) 

DIMENSION  CMEAS (NMIC) ,  MCOM(ll) . 

CMPLEX  X(2-NPFW),  Y(2-NPFW),  W0RK(2-NPFW) ,  XX(2-NPFW),  YY(2-NPFW),  ZZ(2‘NPFW) 

2.  Data  file  Format.  The  input  to  the  program  was  data  taken  from  the 
PASS  Experiment  that  had  been  transferred  to  files  in  the  UNIVAC  1108  system  by  PSL 
personnel.  The  format  modified  slightly  to  include  MET  data  and  file  identi¬ 
fication  for  these  data  files  is  illustrated  in  Figure  2.  Each  data  file  is 
728  lines  in  length  with  the  first  eight  lines  containing  timing  and  shot  in¬ 
formation  while  the  remaining  lines  contain  the  data  for  each  microphone.  To 
be  more  specific  in  the  format  we  have  the  following  line  by  line  structure. 


12  12 

13 

12 

12 

F6.3 

F4.1 

14 

12 

— 

— 

— 

Line  1: 

_ 6 _ 1 

341 

12 

36  _ 

2.024 

_  06 . 7 

C380 

09 

a  b 

c 

d 

e 

f 

g 

h 

i 

a  =  Number  of  microphones 
b  =  PASS  source  number 
c  =  Day  of  Pass  shot 
d  =  Hour  of  Pass  shot 


e  =  Min  of  Pass  shot 
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Line  2 : 


Lines  3-8 


Line  9: 


Line  10: 


Line  128: 


f  =  Second  of  Pass  shot 

g  =  Effective  temperature  in  Deg  C.  for  shot 
h  =  Wind  direction  in  MILS 
i  =  Wind  speed  in  KNOTS 
12 

6  SOURCE  1  MIKE  ARRAY  B _ FILE  (S1-D341-MB-1 . ) 

a  b 

a  =  Number  of  microphones 

b  =  Indication  of  source  position,  microphone  array  and  file  identi- 
f icat ion 

12  14  13  12  13  12  12  F6.3 

: _ 1 _  1920  _ 120  _  27  _  341  _  12  _  36  _  35.611 

a  b  cdefg  h 

a  =  Microphone  index  J 

b  =  Number  of  values  per  microphone  in  record 
c  =  Number  of  lines  of  data  for  each  microphone 
d  =  PASS  microphone  number:  MIC(J) 
e  =  Starting  day  of  microphone  record 
f  =  Starting  hour  of  microphone  record 
g  =  Starting  min  of  microphone  record 
h  =  Starting  second  of  microphone  record 
F5 . 0  F5 . 0  F5 . 0 

_ 12 _ 14  •  •  • _ 15 

al  ai  a16 

=  value  of  microphone  1  at  time  i  =  1,...,  16 
a^  =  value  of  microphone  1  at  time  i  =  17,...,  32 

a.^  =  value  of  microphone  1  at  time  i  =  1905,...,  1920 


In  a  similar  fashion 


129-248: 

values 

of 

microphone 

2 

at 

times 

i  =  1,  2,.., 

.  ,  1920 

249-368: 

values 

of 

microphone 

3 

at 

times 

i  =  1,  2,.., 

.  ,  1920 

369-488: 

values 

of 

microphone 

4 

at 

times 

i  =  1,  2,.., 

.  ,  1920 

489-608: 

values 

of 

microphone 

5 

at 

times 

i  =  1,  2,.., 

.  ,  1920 

609-728: 

values 

of 

microphone 

6 

at 

times 

i  =  1,  2,.. 

.  ,  1920 

3.  Bringing  in  Data  from  Master  File.  In  the  original  version  of  the 
program  a  subroutine  was  provided  by  PSL  to  read  the  data  files  and  transfer 
this  information  to  the  SMIC(J,I)  indexed  array.  In  the  present  version  the 
new  subroutine  RDATA(SMIC, MIC, IDAY , IHOUR, IMIN , SEC, TEMP, MILS, KNOTS)  was  written 
to  obtain  a  desired  output  format  and  include  meteorological  and  other  cata¬ 
loging  information.  A  Fortran  listing  of  the  subroutine  appears  in  Appendix  B. 
To  bring  in  the  data  to  the  program  arrays  simply  requires  the  Fortran  state- 
men  t : 

_  CALL  RDATA(SMIC, MIC, IDAY, IHOUR, IMIN, SEC, TEMP, MILS, KNOTS) 

The  indexed  arrays:  SMIC,  MIC,  IDAY,  IHOUR,  IMIN,  SEC  are  described  in 
the  previous  section  and  represent  the  values  and  starting  time  information 
for  each  microphone  signal,  while  the  TEMP,  MILS,  and  KNOTS  are  defined  below. 
TEMP  =  Effective  temperature  at  the  time  of  the  shot  in  degrees  Centigrade 
MILS  =  The  angle  of  wind  in  mils  at  the  time  of  the  shot 
KNOTS  =  The  speed  of  the  wind  at  the  time  of  the  shot 
B .  Determination  of  Time  Windows  and  Energies 

In  many  cases  the  time  intervals  of  the  received  signals  are  much  larger 
than  the  transient  signal  searched  for  and  if  the  correlation  procedure  is  used 
on  the  entire  length  of  each  signal  an  exorbitant  amount  of  time  would  be  con¬ 
sumed.  Therefore,  windows  are  established  around  the  location  of  the  signal 
bursts.  This  part  of  the  procedure  consists  of  computing  rough  arrival  times, 
starting  points  for  windows  and  the  energies  within  the  windows. 
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1.  Computation  of  Rough  Arrival  Times.  The  rough  arrival  times  are  deter¬ 
mined  by  the  use  of  an  energy  detector  for  each  signal.  That  is  the  energy 
within  a  sliding  time  window  (I,  I+NSW)  is  calculated  for  all  values  of  time  I 
equal  1  to  NL-NSW.  The  point  in  time  1^  for  the  Jth  microphone  signal  where 

the  energy  is  maximum  is  called  the  rough  starting  time  LRST(J).  If  X(J,I) 

I  =  1,2,...,NL  is  the  input  to  the  energy  detector  for  the  Jth  signal,  the 
running  value  of  the  energy  within  the  Jth  sliding  window  becomes 

I+NSW 

E (J ,  I)  =  l  SMIC2 (J ,k)  I  =  1,2 . NL-NSW 

k=I 

LRST(J)  =  value  of  I  such  that  E(J,I)  is  maximized. 

2.  Computation  of  Fourier  Window  Starting  Points.  The  starting  point 
NR(J)  for  the  Jth  signal  Fourier  window  is  selected  by  bracketing  the  rough 
starting  time  on  both  sides  by  NSW/4.  Care  must  be  taken  to  make  sure  that 
the  beginning  and  ending  points  of  the  window  remain  within  1  and  NL.  This 
is  done  by  making  the  following  assignment  for  NR(J). 

1  if  LRST(J)  -  NSW/4  <  1 

NR(J)  =  LRST(J)  -  NSW/4  if  1 <  (LRST(J)-NSW/4)  i  NL  -  3NSW/2 

NL-3-NSW/2  if  (LRST(J) -NSW/4)  >  NL  -  3NSW/2 

3.  Computation  of  Signal  Energy  within  the  Fourier  Window.  Once  NR(J) 
has  been  determined  the  energy  within  each  Fourier  window  can  be  calculated  as 
follows : 

NR(J)+3NSW/2 

E(J)  =  l  SMIC2 (J ,k) 
k=NR(J) 

The  starting  point  of  the  Fourier  window  and  signal  energy  are  then  printed 
out  for  each  microphone. 

C.  Prefiltering  Operation 


In  many  cases  there  is  noise  present  on  the  lines  during  the  duration  of 
the  signals.  The  noise  could  come  from  a  wide  variety  of  sources  including 
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wind  noise,  sixty  cycle  interference,  and  machine  related  noise.  To  charac¬ 
terize  all  these  noises  analytically  becomes  almost  impossible.  The  program 
was  designed,  however,  to  try  to  minimize  the  effects  these  noises  would  have 
by  prefiltering  the  data.  This  was  done  by  using  estimates  of  the  spectrum 
of  the  noise,  and  the  spectrum  of  the  signal  plus  noise. 

1.  Determination  of  Starting  Point  for  Noise  Window.  The  noise  window 
for  each  signal,  selected  to  be  before  the  Fourier  window  which  contains  the 
signal  plus  noise,  consists  of  512  points  consistent  with  the  Fourier  window 
in  size.  By  keeping  the  windows  the  same  size  the  determination  of  the  pre¬ 
filter  and  application  of  the  prefilter  is  simplified.  The  starting  point  for 
the  noise  window  for  each  signal  was  determined  by 

Start  Point  (J)  =  NR(J)  -  512 

If  this  was  less  than  zero  for  some  J  then  the  starting  point  for  the  Jth  noise 
window  was  set  equal  to  the  first  data  point  of  that  signal  record.  Since  most 
of  the  data  records  had  over  500  samples,  i.e.  500  msec  worth  of  noise  preced¬ 
ing  the  signal,  few  cases  of  signal  and  noise  window  overlaps  were  reported. 

2.  Determination  of  Estimate  of  Noise  Spectrum.  For  each  pair  of  sig¬ 
nals  SMIC(I,L)  and  SMIC(J,L),  and  respective  noise  window  starting  points  KXX 
and  KYY,  a  rough  estimate  of  the  noise  spectrum  is  desired.  The  estimation 
procedure  used  was  that  of  obtaining  the  spectrum  for  each  signal  and  simply 
averaging  the  results.  To  provide  a  spectrum  compatible  in  samples  the  FOURG 
SUBROUTINE  was  used  on  each  of  the  following  signal  vectors. 

512  samples  512  samples 

/■ - * - ^  ✓ - A - V 

XX  =  [SMIC(I ,KXX) ,  SMIC(I ,KXX+1) , . . . ,  SMIC(I ,KXX+511) ,  0,  0,...,  0] 

YY  =  [SMIC(J.KYY) ,  SMIC (J ,KYY+1) , . . . ,  SMIC(J .KYY+511) ,  0,  0 . 0] 

The  addition  of  zeros  will  not  alter  the  spectrum  determined  and  provide  the 
proper  frequency  indices.  The  spectrum  was  then  estimated  by 
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„  ,  .  FOURC(XX)  •  FOURG* (XX)  +  FOURG(YY)  •  FOURC* (YY ) 

Spectrum  for  noise  = - 1 - - - ^ - 1 — — - - — - 

The  noise  spectrum  was  stored  back  in  the  real  part  of  the  complex  array 
YY(L) ,  L  =  1,2,...,  1024. 

3.  Determination  of  Estimate  of  Signal  Plus  Noise  Spectrum.  For  each 

pair  of  signals  SMIC(I,L)  and  SMIC(J,L)  and  respective  Fourier  window  starting 

points  NR(I)  and  NR(J)  the  following  signal  vectors  were  defined. 

512  samples  512  samples 

- - * - ^  , - « - , 

XX  =  [ SMIC (I ,NR(I) ) ,  SMIC (I ,NR(I )+l) ,. . . ,  SMIC(I ,NR(I)+511) ,  0,  0,...,  0] 

ZZ  =  [SMIC(J ,NR(J) ) ,  SMIC(J,NR(J)+1) ,...,  SMIC(J ,NR(J)+511) ,  0,  0,...,  0] 

The  signal  plus  noise  spectrum  was  stored  back  in  the  real  part  of  complex 

array  XX (L) ,  L  =  1,2,...,  1024. 

4.  Computation  of  Prefilter  Values.  The  prefilter  that  is  to  be  applied 

to  each  pair  of  signals  was  selected  to  be  of  the  form 

H(L)  =  (L)/(4  (L)  +  *  (L))  L  =  1,2,...,  1024 

ss  no 

where  (L)  and  represent  Che  value  of  the  energy  spectrum  density  of 

signal  and  noise  respectively  at  index  L.  Since  the  0gs(L)  was  not  known  the 
prefilter  frequency  response  was  written  in  the  following  equivalent  form 

H(L)  =  1  -  0nn(L)/($ss(L)  +  0nn(L» 

Since  only  the  estimates  of  $  (L)  and  <5>  (L)  +  <5>  (1)  are  available,  care 

nn  ss  nn 

must  be  taken  to  insure  that  H(L)  does  not  go  negative.  The  program  uses 

Real  (XX(L))  and  Real  (YY(L)>  to  represent  0  (L)  +  4>  (I.)  and  $  (L)  re- 

ss  nn  nn 

spectively  and  defines  PRC  and  PRD  as 

PRC  =  REAL(YY (L) ) /REAL (XX (L) ) 

PRD  =  REAL (XX (L)) 

The  prefilter  value  ZZ(L)  at  frequency  index  L  is  described  as  follows  for 
L  =  1,2,...,  1024. 
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ZZ(L)  = 


(0,0) 

(1,1) 

CMPLX(H(L) ,H (L) ) 


If  PRC  >  1 
If  PRD  <  .001 
otherwise 


5.  Application  of  Prefilter.  Once  the  prefilter  values  have  been  deter¬ 
mined  and  the  product  of  the  DFT  of  each  positioned  signal  is  eventually  put 
in  Y(L)  L  =  1,2,..., 1024  the  prefiltering  operation  is  done  simply  by  multi¬ 
plying  each  frequency  component  by  ZZ(L)  that  is 

ZZ (L)  =  CMPLX  (PREF.PREF) 

Y(L)  =  Y(L)*ZZ(L) 


6.  Application  of  Comb  Filter  for  60  HZ  and  Harmonics 
The  power  spectral  density  was  obtained  for  microphone  noise  signals 
from  the  PASS  experiment  showing  dc  power,  power  in  60  HZ  and  harmonics,  some 
low  frequency  wind  noise,  and  some  unidentified  interference  of  approximately 
7.8  HZ  and  higher  harmonics.  An  example  appears  in  Appendix  C.  To  minimize 
the  effects  of  these  interferences  a  comb  filter  including  a  low  pass  filter 
was  applied  in  the  frequency  domain.  For  the  1024  point  Fourier  transform  the 
frequency  indexes  and  corresponding  analog  and  digital  frequencies  are  given 
by : 

index  digital  frequency  corresponding  analog  frequency 

k  (k-1)  •2tt/1024  radians  (k-1)  •  2tt/1024T  radian/sec 

Since  the  time  T  between  samples  is  1  msec  for  the  PASS  data  the  corresponding 
analog  frequencies  are  spaced  2tt/1.024  rad/sec  or  .97656  HZ  apart.  The  D.C. 
component  index  is  one.  Since  60  HZ  and  harmonics  do  not  fall  exactly  on  in¬ 
dex  numbers  this  causes  a  leakage  to  a  band  of  frequencies  that  must  be  con¬ 
sidered.  The  apparent  index  given  in  Table  1  is  determined  by: 

Apparent  index  =  freq. (1.024)  +  1. 

The  uoper  index  band  required  because  of  symmetry  of  the  DFT  is  determined  by 


ku  =  2*NPFW  +  2  -  k  =  1026  -  k 


J. 
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TABLE  1 

Comb  Filter 

Frequency  Bands 

Identification 

Frequency 

Band 

Apparent 

Index 

Lower  Index 
Band  k 

Upper  Inde: 
Band  k 

u 

a 

D.C. 

1 

1 

1 

b 

0+-7.8  HZ 

1-8.9 

2-10 

1016-1024 

c 

60  HZ 

62.44 

60-65 

961-966 

d 

120  HZ 

123.88 

121-126 

900-905 

e 

180  HZ 

185.32 

183-188 

838-843 

f 

140  HZ 

246.76 

244-249 

777-782 

The  Index  information  for  the  comb  filter  was  given  in  the  following  data  card 
_ DATA  MCOMB/1,  2,  10,  60,  65,  122,  126,  183,  188,  244,  249/ 

V— ’ 

f 


e 


ah  c  d 

The  comb  filter  was  applied  by  setting  the  values  corresponding  to  the  index 
bands  equal  to  zero  in  the  frequency  domain  output  of  the  pair  by  pair  corre 
lation. 

D.  Pair  by  Pair  Correlation 

The  pair  by  pair  cross  correlation  R  (L)  of  the  windowed  signals  is 

xy 

accomplished  in  the  frequency  domain  by  use  of  the  FFT  algorithm.  That  is 


R  (L)  =  IDFT[DFT(X(L)) -DFT*(Y(L))  ] 
xy 

Since  the  above  operation  performs  circular  correlation  it  is  necessary  to 
position  the  data  signals  X(L)  and  Y(L)  such  that  the  circular  and  linear 
correlation  are  the  same.  The  *  means  complex  conjugate  and  is  needed  to 
provide  correlation  rather  than  convolution. 

1.  Positioning  of  Signals  for  Determining  the  Cross  Correlation.  To 
obtain  linear  convolution  resulting  in  proper  correlation,  the  signals  SMIC(I,L) 
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and  SMIC(J,L)  L  =  1,2,..., 512  are  positioned  as  follows 

I,  1  ••  2  •••  512  513  514  1024 

Real  X(L)  SMIC(I ,1) ,SMIC(I ,2) , . . . ,SMIC(1 ,512) ,  0  0  •••  0 

I mag  X(L)  0  0  ...  0  0  0  •••  0 

Real  Y(L)  0  0  ...  0  SMIC (J , 1) , SMIC ( J , 2) , . . . , SMIC (J , 512) 

Imag  Y (L)  0  0  ...  0  0  0  •  •  •  0 

2.  Finding  the  Raw  Time  Differences  between  Windowed  Signals.  Once  the 

DFT(x (L ) ) • DFT*(Y (L) )  is  prefiltered  and  comb  filtered  the  cross  correlation 

R  (L)  is  obtained  back  in  the  real  part  of  array  Y (L)  L=1 ,2 ,  . . . ,1024  by  taking 
xy 

the  Inverse  Fourier  Transform.  The  Fortran  statement  uses  a  +1  to  indicate  in¬ 
verse  as  follows. 

_  CALL  FOURG (Y , 1024 ,+l ,  WORK) 

A  search  is  then  performed  on  | Real  (Y(L) ) | ,  L=l,2 . 1024  to  obtain  the  maximum 

value  CMAX  of  the  cross  correlation  and  the  time  TMAX  of  the  maximum  value.  Be¬ 
cause  of  the  positioning  of  the  signals  the  raw  time  difference  is  then  given  by 

TAU(I , J)  =  TMAX  -  NPFW  -  1 

3.  Determination  of  the  Normalized  Correlation  Coefficient.  The  normalized 
correlation  coefficient  GAM(I,J)  is  obtained  by  dividing  the  maximum  value  by  the 
square  root  of  the  product  of  the  energies  and  is  given  by 

GAM(I , J)  =  CMAX/((E(J)*E(I))  **0.5) 

If  the  signals  [SMIC(I,L)]  and  [SMIC(J,L)]  are  identical  in  shape  but  of  dif¬ 
ferent  amplitudes  the  GAM(I,J)  will  be  equal  to  1.  If  there  shapes  are  signi¬ 
ficantly  different  the  value  of  GAM(I,J)  will  be  approximately  zero.  In  this 
way  the  normalized  correlation  coefficient  provides  a  measure  of  similarity  of 
the  two  signals  tested. 

4.  Determination  of  the  Rough  Time  Differences  between  Signals  I  and  L. 

The  raw  time  difference  represents  the  time  difference  between  the  signals  in 
the  windows.  To  obtain  the  rough  time  differences  between  microphone  signals 
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the  raw  time  difference  TAL’(1,1.)  must  be  adjusted  by  the  starting  times  of  the 
windows  as  follows 

l\TAl'(  1,1.)  =  FLOAT (IT(I)-IT(L)  ]  +TAU(I,L) 

These  roue.lt  tint.'  differences  TAU(I,L)  1  =  1 , 2  ,  .  .  .  ,NMIC  ,  and  L=I ,  1+1  ,  .  ,  .  ,NMJC 
are  part  ot  the  in,.ut  to  the  least  square  adjustment  section. 

L .  Determination  of  Least  Square  Time  Fit 

The  rough  time  difference  RTAU(I,J)  and  normalized  correlation  coefficient 
GAM(I,.!>  1  =  1 ,2  , . . .  ,NMI<’  J=1 ,2 , . . . ,NMIC  are  now  used  to  obtain  a  consistent  set 
of  time  differences  by  using  a  weighted  least  squares  procedure.  Each  step  in 
the  procedure  is  described  in  the  following  paragraphs. 

1.  Establishment  of  Weights  for  Least  Square  Time  Fit.  A  weight  is  at¬ 
tached  to  each  one  of  the  rough  time  differences  indicating  a  measure  of  relia¬ 
bility  of  those  estimates.  The  measure  or  weight  G(T,J)  assigned  to  RTAU(I,J) 
is  defined  to  be  a  function  of  the  normalized  correlation  coefficient  GAM(I,J). 
Different  functional  relationships  have  been  played  with,  but  at  present  no  way 
of  favoring  one  over  the  others  has  been  established.  The  program  in  the  present 
form  uses  the  f.  flowing  weighting. 

G(I,J)  =  (GAM(I , J) ) 2  if  GAM ( I , J )  <  0.5 

=  ( GAM ( I , J ) ) 1 ^ 2  if  GAM(I , J)  >  0.5 

This  attaches  slightly  more  influence  to  time  difference  values  that  have  a 
normalized  correlation  coefficient  greater  than  0.5  over  those  that  have  cor¬ 
relations  less  than  0.5. 

2.  Determination  of  Least  Square  Time  Fit  with  MIC  1  as  a  Reference. 

The  weights  G(I,J)  described  above  are  used  to  define  a  measure  of  performance 
e  given  by 

NMIC  NMIC 

e  =  l  l  G(I , J)  (RTAU (I , J)  -  (TIME(I)  -  TIME(J)))2 
I>J  J=1 
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The  objective  is  to  find  the  times  TIME(l)  1  =  1,2 . NMIC  such  that  e  is  mini¬ 

mized.  A  solution  of  this  problem  for  G(I,J)=1  is  given  by  the  Author  [3]  and 
requires  solving  for  n  microphone  signals  a  set  of  simultaneous  linear  equations 
in  n-1  unknowns  when  a  given  microphone  has  been  selected  as  a  reference.  in 
Appendix  C  a  solution  for  arbitrary  G(1,J)  with  MIC1  as  a  reference  is  presented 
resulting  in  the  following  set  of  simultaneous  equations  in  t^  where  t^  is  the 
relative  time  of  arrival  of  the  ith  signal. 


G  (1 , 1)  G(2 , 3)  G(2 ,4)  G(2 ,5)  G(2,6)  t> 

G (2 , 3)  G (2 , 2)  G ( 3 , 4 )  G(3,5)  G(3,6)  t2 

G(2 ,4)  G(3,4)  G(3,3)  G(4,5)  G(4,6)  t3 

C.  ( 2 , 5 )  G  ( 3 , 5 )  G(4 ,5)  G(4,4)  G(5,6)  t4 

G(2,6)  G(3,6)  G(4,6)  G(5,6)  G(6,6)  t5 

where 

G(1 , 1>=-(G(1 ,2)+G(2 , 3)+C(2  ,4)+G(2  ,5)+G(2 ,6)) 

G(2 ,2)=-(G(l ,  3)+G(2 , 3)+G(3,4)+G(3 ,5)+G(3,6) ) 
C.(3,3)=-(G(1,4)+G(2,4)+G(3,4)+G(4,5)+G(4,6)) 
G(4,4)=-(G(1,5)+G(2,5)+G(3,5)+G(4,5)+G(5,6)) 
G(5,5)=-(G(1,6)+G(2,6)+G(3,6)+G(4,6)+G(5,6)> 

and 

a2  =  -  C(1,2)R(1,2)+G(2,3)R(2,3)+G(2,4)R(2,4)+G(2,5)R(2,5)+G(2,6)R(2,6) 

a3  =  -  G(1,3)R(1,3)-G(2,3)R(2,3)+G(3,4)R(3,4)+G(3,5)R(3,5)+G(3,6)R(3,6) 

a4  =  -  G(1,4)R(1,4)-G(2,4)R(2,4)-G(3,4)R(3,4)+G(4,5)R(4,5)+G(4,6)R(4,6) 

as  =  -  G(1 ,5)R(1 ,5)-G(2 ,5)R(2 ,5)-G(3 ,5)R(3,5)-G(4 ,5)R(4 ,5)+G(5 ,6)R(5 , 6) 

a  -  -  G(1,6)R(1,6)-G(2,6)R(2,6)-G(3,6)R(3,6)-G(4 ,6)R(4 ,6)-G(5 ,6) R(5 ,6) 

6 

with 


R(I , J)  =  RTAU(I , J) 
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This  set  of  equations  was  solved  by  using  the  CJR  subroutine  package  from  the 
1108  Large  Scale  Systems  Math  Pack  which  is  fully  described  in  | 5 ] .  The  sub¬ 
routine  is  accessed  with  the  following  call 

_  CALL  CJR(A ,6 , 5 , 5 ,6 , $400 , JC , V) 
abcde  f  gh 

where 

a  =  augmented  coefficient  matrix 
b  =  maximum  number  of  columns  in  A 
c  =  maximum  number  of  rows  in  A 
d  =  the  number  of  rows  in  A 
e  =  the  number  of  columns  in  A 

f  =  statement  number  control  is  returned  if  an  overflow  is  detected 
g  =  control  array 

h  =  on  input  V(l)  is  the  option  indicator,  which  to  solve  equations  is  4, 
on  output  V(2)  contains  the  value  of  the  natural  log  of  the  absolute 
value  of  the  determinant  and  V(l)  contains  the  sign  of  the  determi¬ 
nant 

On  output  the  last  column  of  A  is  the  solution  vector  [ 1 2 , 1 3 . . .  t  ]  and  is 
called  AT(2) ,AT(3) , . . . ,AT(6),  Since  microphone  1  was  defined  to  be  the  reference, 
AT(1)  is  set  equal  to  zero.  The  AT[I]  represents  the  starting  time  for  the 
signal  in  the  Ith  window  relative  to  the  signal  in  window  1. 

3.  Adjustment  of  Time  Differences.  Using  the  AT  vector,  a  raw  corrected 
time  difference  can  be  found  by  the  difference  of  the  components 

TAU(I , J)  =  AT(I)  -  AT(J)  1  =  1,2,.. .,6  I<  J<ft 

The  least  square  corrected  time  difference  in  milliseconds  can  then  be  formed  by 
CTAU(I , J)  =  TAU( I , J)  +  FLOAT [IT (I) -IT (J) ] 

4.  Selection  of  MIC  to  be  used  as  reference.  To  select  a  microphone  as 

a  reference  for  determining  the  relative  times  for  the  position  estimator  requires 


a  specification  of  a  performance  criterion.  We  would  like  to  select  the  MIC 
as  reference  that  gives  the  best  possible  relative  times  for  the  position  esti¬ 
mator.  Upon  examining  data,  it  was  found  that  position  estimates  were  sensi¬ 
tive  to  the  selection  of  a  MIC  reference  even  though  a  least  squares  precedure 
had  been  performed  on  the  time  differences.  No  direct  relationship  was  es¬ 
tablished  however  and  the  following  procedure  represents  a  reasonable  wav  to 
select  the  reference  with  no  sense  of  optimality  implied.  As  i  he  normalized 
correlation  coefficient  is  a  measure  of  reliability  of  the  time  differences 
determined,  the  sum  of  the  correlation  coefficients  for  all  time  differences 
associated  with  a  particular  microphone  represents  a  measure  of  reliability 
of  the  time  differences  for  that  microphone.  This  measure  was  called  CMEAS(I) 
1=1, 2,..., 6  and  given  by 

6 

CMEAS(I)  =  y  0(1, .1) 

J  =  1 

The  value  K  for  which  CMEAS(K)  K=l,2,...,6  is  a  maximum  is  called  KMAX  and 
represents  the  number  of  the  microphone  to  be  used  as  a  reference.  In  the  out¬ 
put  of  the  times  for  the  position  estimator  the  reference  microphone  is  identi¬ 
fied  by  having  the  time  0.000. 

5.  Determination  of  Relative  Times  for  the  Position  Estimator.  The 
position  estimator  USRAN1[1]  requires  six  relative  arrival  times  and  the 
proper  meteorological  data  as  input  to  determine  the  position  of  the  source. 

The  number  (KMAX)  of  the  microphone  signal  that  was  selected  as  a  reference 
was  determined  previously.  To  find  the  times  relative  to  this  microphone 
could  be  easily  accomplished  by  filling  out  the  lower  triangular  part  of  the 
CTAU(I,L)  array  with  the  negative  of  upper  triangular  part.  In  this  way  the 
time  differences  with  respect  to  each  microphone  appear  on  the  rows  of  CTAU(I,L) 
and  the  selection  of  the  KMAX  microphone  corresponds  to  selecting  a  row. 
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F .  Printing  the  Output 

An  example  of  the  overall  program  output  format  is  given  in  Figure  j  and 
consists  of  the  four  basic  parts  described  below. 

1.  Shot  Information.  The  header  information  is  in  reference  to  the  data 
files  from  the  PASS  experiment  and  gives  the  day,  time,  and  position  of  the 
shot  as  well  as  the  position  number,  starting  times  and  lengths  (MSEC)  of  the 
data  record  for  the  microphones  recording. 

2.  Window  Information.  The  window  starting  time  in  MSEC  determined  by 
the  program  relative  to  the  starting  times  for  each  microphone  record  are 
given  along  with  the  signal  energy  within  the  window.  The  microphones  numbered 
1  through  6  correspond  to  the  normal  ordering  of  the  PASS  microphone  numbers. 
The  window  signal  energy  and  starting  time  could  be  used  in  an  interactive 
mode  (not  programmed  at  this  time)  to  allow  the  operator  to  check  if  a  signal 
appears  in  the  window  or  not. 

3.  Time  Differences  Information.  The  normalized  correlation  coefficient, 
time  differences  in  MSEC,  and  corrected  time  differences  in  seconds  are  given 
for  each  pair  of  microphone  signals.  These  results  may  be  used  to  infer  the 
overall  reliability  provided  by  the  normalized  correlation  coefficient  and  a 
measure  of  consistency  by  examining  the  changes  made  in  the  rough  time  dif¬ 
ferences  RTAU(I,J)  to  get  the  corrected  time  differences. 

4.  Met  and  Timing  for  Position  Estimator.  The  times  in  seconds  and  the 
Met  information  that  could  be  used  as  input  to  a  position  estimator  are  pro¬ 
vided  for  interfacing  purposes.  The  timing  information  appears  in  an  array 
TIMF.(J)  J=l,2,...,6  while  the  temperature,  wind  direction  and  wind  speed  are 
labeled  TEMP,  MILS,  and  KNOTS  respectively. 

III.  USING  THE  PROGRAM 

The  procedure  given  in  this  section  will  apply  to  the  use  of  the  program 
on  the  WSMR  UNIVAC  1108  system.  Presently  the  program  is  stored  under  file 
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M  LC ( 32  >  I 

h: IS  HRS 

39 . 4S  l 

sr  c 

rut-  WINUUW  STARTING  TIMES* MSEC >  RELATTVt  TO  APOVE 

MIC  l  MIC  2  MIC  .5  MIC  4  MIC  '5  MIC  6 

742  7  26  A/7  A  6  6  686  5  A 5 


me  SICNA I.  CNII-Gi  W.L  fl-fJiV  EACH  UfNDriU 


MIC  1 

m  i  c  2 

MIC  3 

MIC  A 

Mic  ■: 

3  31,9. 

264966  . 

H322272. 

21 I 12 423. 

3541 5047. 

34 

:0I; l\!i !  AT  TON. 

r  t  me;  nil  r 

“LRLNCES ’ MSEC  > 

*  AND  CORIM 

ZCTEO  TIME  RIFE 

ERENCC.: 

!t£N  EACH 

PAIR 

OF  M I CROPHONF  ->  1 1  •I'MUi 

GAMMA < 

!. .  2  )  - 

.693 

RT  AIM  l  *  ?.  ' 

5 . 0 

C  T  A  U  (  1  *  2  >  “ 

.014 

GAMMA 

!  *  > ) 

.  703 

R T  nU (1*3) 

-4  22.0 

CTAU< 1  *  5 ‘  - 

-.436 

GAMMA  * 

1  rl  1  - 

,  6o  V 

RTAtK  .1*4) 

*-1321.0 

CTAU (1*4) 

1.340 

GAMMA < 

!  *  5  i  - 

.  ASS 

RTAU( 1 *5) 

--2 AS 2.0 

C  TAU<  1  •  5  )  -- 

-2.A42 

GAMMA ‘ 

1  •  O  >  := 

.  7  AS 

RfAU<  1.  *6> 

=  4329.0 

CTAU  <  1  *o)  ■- 

-4.319 

GAMMA  •. 

•  t  \  )  ~ 

.A  74 

RT  AU ( 2 , 3  > 

=  -455.0 

CTAU (2*3)  = 

-.451 

Gil  MM  A  ( 

::  -  4  >  - 

.703 

R.TAI.K2*  4  ) 

--  1359.0 

C  T  AU (2*4) - 

-1.354 

GAMMA * 

J  r  5  > 

.  730 

F:  TAi  «'2*S) 

- -CASA « 0 

C  T  AU  (  2  *  5  )  ■= 

-2.656 

GAMMA < 

2  *  •',  )  ‘ 

.753 

R  TAU <2*6) 

••••-4353 . 0 

CTAU  (  2  *  6  >  ■= 

-4.333 

G  A  Mil  A  * 

5*4)  ••• 

.7  74 

R  TAU  <  3 .  4  ) 

-  -  90  3.0 

CTAU (3*4) - 

-.903 

GAMMA • 

S-S )  - 

♦  7?>6 

RT  Al)  (3*5) 

--2201 .0 

CTAU (3*5) = 

-2.206 

GAMMA  < 

3  *  A  )  ■■■ 

.72  4 

R TAIJ(  3-6) 

=•-3377 . 0 

CTAU ( 3  *  A )  = 

-3.882 

GAMMA  ■’ 

4*5  >  -- 

♦  ’95 

RTAU  \  4*5) 

•=-  1298.0 

CTAU (4 *5) - 

-1.303 

GAMMA  < 

4  *  6  )  - 

.  75  5 

R  r  AU < 4  *  A ) 

*-2974.0 

C  T  AU (4*6)  - 

-2.9  7? 

GAMMA  < 

S*A) 

.734 

RTAU (5* 6> 

= - 1 6  7  A  .  0 

CTAU (3 *6) = 

-i .676 

rill:  I  IMIS  IN  SEC  rn  £<£  USEtl  FOR  INPUT  TO  A  POSIT  lOil  estimator 


TIME! J. )  =• 

T  IMG  ( 2 )  -- 

t  r  me ( i > - 
r  J HE (  4  )  - 
TliiE'O)  = 
TIME (A  )  = 


-4.319  SEC 
-4.3  13  SEC 
-3.082  SEC 
-2.979  SEC 
-  I  . A  7  A  SEC 
.000  SEC 


rill  Ml.  r  ro  i/|  ;i  ,|  (I  INPUT  TO  POSITION  ESTIMATOR 

TEMP  5.4  DC03  C. 
WIND  DIM  C T1UN =  A1S0  MILS 
WIND  SPEED  -  9  KNOTS 


Figure  3 


Example  of  the  program  output  format 
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name  and  element  LCLPF.SOUNDR  in  both  a  symbolic  and  absolute  form.  The 
following  control  statements  will  execute  the  program; 

@ASC,  AZ  S3-D341-MB-1 .  Assign  data  file  to  run 

@USE  12.,  S3-D341-MB-1  Use  file  S3-D341-MB-1  as  logical  unit  12 

@XQT  LCLPF.SOUNDR  Execute  program. 

If  the  program  is  to  be  executed  on  another  computer  a  duplicate  of  the 
symbolic  form  of  SOUNDR  must  be  obtained  and  if  a  different  data  structure  is 
to  be  used  the  subroutine  RDATA  must  be  rewritten.  The  main  program  and  sub¬ 
routine  can  then  be  compiled  on  the  computer  to  be  used. 

IV.  RESULTS 

The  program  LCLPF.SOUNDR  was  run  using  various  signals  from  the  PASS  ex¬ 
periment  as  input.  The  results  for  the  test  shots  given  in  Table  2  are  shown 
in  Figures  4  through  10.  The  times  and  Met  data  indicated  in  Figures  4  through 

10  were  used  in  a  position  estimator  program  resulting  in  the  Miss  distances  given 
in  Table  3. 

TABLE  2.  PaSS  Data  Shots  Used  in  Program  Testing 
UNIVAC  1108  Identification 


File  No. 

Source 

Day 

MIC  Array 

Number 

S1-D341-MB-1. 

1 

341 

B 

1 

S1-D341-MB-2 . 

1 

341 

B 

2 

S2-D341-MB-1. 

2 

341 

B 

1 

S2-D341-MB-2. 

2 

341 

B 

2 

S2-D341-MB-4. 

2 

341 

B 

4 

S3-D341-MB-2. 

3 

341 

B 

2 

S7-D341-MB-2. 


7 


341 


B 


2 
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day  in:  shot  rihi:  'Ji.u,  hrs  sue. 

Sl'IJIVCE  l  rllKL  array  n 

mii:s:  77.23. 29  . io. ,;i . 3::  i  v ; “ <>  vauils/mic 

i m  siautiNu  i  (mi  s  r ui;  each  m [cmipiiunl  gicnai 


nil  \2/> 

:  r i :  36  mi <8 

35  .  A 1 1 

SEC 

M  2  C ( 30  )  : 

r. :  36 

HRS 

39. 101 

SEC 

m  n:<28) 

:  5J36  hks 

36 . SO 7 

SEC 

MIC (31 > : 

v» :  36 

HRS 

11 . 199 

SI  r 

MIC (29) 

:  5 : 16  hrs 

37*78/ 

811 C 

MIC (32) : 

5 :  36 

HRS 

1 3 . 803 

SEC 

WINDOW 

STARTING  T I MI¬ 

S ( MSEC ) 

kE.I 

-ativc 

to  abooi: 

MIC  1 

MIC  2 

MIC  3 

MIC  4 

MIC  *3 

M I L  6 

808 

603 

6  GO 

71  A 

701 

7  A  7 

SIGNAL  Li 

'i-RGY  WITHIN  EACH 

WINDOW 

MIC  1 

MIC  2 

MIC  3 

MIC  1 

3622710. 

->073-12  t 

9752818. 

26270111 

MIC  5  MIC  6 

I140V67C.  1395564.?. 


THE  CORRELATION.  I IMC  DIFFE  RF  NOUS  <  MSEC  >  .  AND  (  ORRECrED  TIME  I']  FT  t  RT  NC.EC  (  SEC  > 
BETWEEN  EACH  PAIR  Of  MICROPHONE  SIGNALS 


GAMMA < it 2>- 

*  553 

R  T AU ( 1 *2> 

=  -884.0 

CTAlKlt2>«  .884 

GAMMA (1,3)  - 

R TAU ( 1 p  3 ) 

“-210 L  .0 

C  TALK  1.3) --2. 182 

GAMMA  (  1  p  4  >  := 

.4/1 

R  TAU (  1 p  4  ) 

-- 3860.0 

CTAUCl  .4)  --3.86J 

GAMMA  ( 1  p  5  > :- 

.11 B 

R I AU ( 1 . 5 ) 

“-GO 7 8.0 

C  TALK  1 .5 >  =  -5.883 

GAMMA  < 1 p  6  )  = 

.534 

In  TAU  ( 3  . 6  ) 

- -021 8.0 

CTAUU  p6)~  —  8.216 

GAMMA  <  2  r  3  >  == 

.  395 

E \ AU (2*3) 

1295.0 

C  T  AU  (  2  p  3  )  -  -  1  .797 

GAMMA ( 2 p  4  >  ™ 

.412 

IV  TAU  (2*4) 

-  -2777,0 

CT AU  <  2  p  4 ) r -  2 . 9 7 7 

GAMMA  <  2  p  5 ) » 

.336 

RT  AU (2*5) 

•--4993.0 

C  T  AU  <  2  p  5 )  ■-  -  4  .  v9 8 

GAMMA <  2  p6  > - 

♦  425 

IV !  AU  ( 2  p  6  ) 

~~ 7334 . 0 

CTAU(2p6)  --7.331 

GAMMA  (3.4)  = 

.421 

RT AU (3*4) 

.-•-1677.0 

CT  A  LI  (3*4)  1  .679 

GAMMA  ( 3  9  5  > :- 

.  370 

R  T AU < 3  *  5 ) 

--3695.0 

C T AU <  3 p 5  > • -3.701 

GAMMA ( 3  p  6 ) - 

.  457 

R TAU (3.6) 

*-6037.0 

CT AU <  3  p  6 ) ~~6 . 034 

GAMMA ( A r 5 ) = 

.422 

R  TAIJ  (4.5) 

-■--2017.0 

CTAU(4pS)- -2.022 

GAMMA (4 » 6 )- 

.439 

R 1 AU ( 4 . 6 ) 

--4356.0 

C  TALK  4*  6)-=  -4 . 355 

GAMMA  <  5  p  6  )  = 

.376 

R  T All  ( • j  p  6  ) 

-  -2310 . 0 

CTAUCSp  6)^  -2 .337 

TIMES  IN  SEC  TO  BE  US EH  FOR  INPUT 

TO  A  POSITION  ESTIMATOR 

T I  ME  ( 1  )  := 

.OOO 

SEC 

T IME  <  2 )  = 

.  884 

SEC 

TIME (3)” 

2.182 

SEC 

T  [  ME  ( 4  )  = 

3.861 

SEC 

TIME <5>“ 

5.833 

SEC 

TIME ( 6 ) ~ 

8.216 

SEC 

THU  MET  TO  A  BE  USED  AS  INPUT  TO  POSITION  ESTIMATOR 

TEMP=  6.7  HCGS  C. 
WIND  DIRECTION-  6380  MILS 
WIND  SPEED-  9  KNOTS 


Figure  4.  Results  from  PASS  shot  S1-D341-MB-1 
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HAY  34  L  :  '.Mil  1  I  1(111  CIIjI.  111  ;..  1.640  SEC 

'(ll  IKCC  I  liIKt:  AI-RnT  [l 

MIC;;:  27. 2IJ.. ’9.30. 31 .32  1  Y20  Vf.UICS/rtlC 

I  me  STARTING  r  I  fit  I  Tl|:  I  .VIII  IV.P.kUITIUNP  1.10NAI 


nic< sn : 
die «.  jh  ) : 
mil  ( :>9) : 

St  53  HRS 
a: vi  hrs 
o:sr»  HU’S 

3b. 22 7  ’ 
36. 123  i 
3/ .40  x  • 

)  r 

re 
)  r: 

mjc:(30>:  n 
MU'  71  )  :  i) 
m  j  i;  ( 32  > :  u 

V;  hks  3v  r 
so  nee  4i,  lie 
\j[,  HRS  13.41V 

Cl  c 

'■.1  c 
Cl  c 

UlriDUU  S 

ART  INS  MM! 

*,  \  MSEC ) 

RL'LA  I  1  VC 

TO  ARMVE 

MIC  1 

m  re  ? 

MIC:  3 

tt  1.  C  4 

MIC  '6 

MIC  6 

8?  1 

797 

360 

872 

894 

oco 

11'!  SIGNAL  ENERGY  UlTI-UN  LACH  WINDOW 


MIC  1 

MIC  2 

MIC  3 

NIC  4 

NIC  5 

MIC  6 

208123. 

49282 . 

1073772 . 

15422260. 

3029/288. 

lCCl*. >8  63 

Till  T.TltVRI.  I.A  I  ION.  TIT1F  DIPT  fREWCrs  '  MSCC )  t  ANIi  CORRECTED  T INC.  D I FFCRE NCES ( SEC > 
Bt.iWEEN  LACH  PAIR  OF  1 1 T r  l-'jl  HONIL  SIGNALS 


b  A  rid  A  <  1  f  2 )  = 

0  <Kj 

(jiil  id  A  (  1  f  3  )  •- 

.  1 30 

GAddA  < 1 • 4 )  - 

.  233 

GAMMA  <  1  >  v> )  -■ 

.206 

GAMMA ( 1 r 6 > * 

.  197 

GAMMA <  2 f 3 > « 

.  160 

GAMMA < 2  9  A ) * 

.  307 

GAMMA < 2  f  5 ) « 

.306 

GAMMA  <  2  r  6 )  ® 

.271 

GAMMA <3 9 4)- 

.  435 

GAMMA <  3 ? 5 > « 

♦  496 

GAMMA <  3  9  6  )  - 

.303 

GAMMAS  4  »5> - 

.933 

GAMMA ( 4  9  6  >  = 

.7  71 

GAMMA  (5y 6 )~ 

.807 

RTAIK1  <2>  ■■■  725.0 

R TAUT  If  3)  -  7259.0 
FT  AIK  1*4)  -391  A. 0 
RTAIK  1  -59 .'4.0 

RTlUKl  79o0.0 

RTP.U<2.3>  = -1308.0 
R  TALK 2. 4 ) =-3034 • 0 
RTAU(2»5> — 0042. 0 
RTAUT2.6)-  7313.0 
RTA1K  3.4) - -1733.0 
RTAU < 3 »3>=- 3740.0 
K’TAUT  3*6  )  -  -6003.0 
RTAU (4. 5 >--2007.0 
RTA0<4 *6 >=-4302.0 
RTAU <5. 6 >=-2295.0 


CTAU(1.2>=  -.321 
CTAUi  I  *3  ■  2. 137 
CTAUU  ,4> -—3.853 
C.TAU(1*5>  —  5.040 
CTAUI  1  *6)  -■  -8.  I  43 
CIAU<2» 3>=-l .316 
t:TA1K  2  *  4  )  -  -3 . 033 
C  T AU ( 2 * 5) -- 5.040 
CTAIK2.6  >--7.327 
CT  AU ( 3  *  4 '  - -  1 . 716 
CTAU(3*5>-  3.723 
C I AU  Co  *  6 ) --6 .011 
CTAUI 4. 5)  --2. 00 7 
CTAU< 4.6) --4.294 
CTAUI 5. 6 >*-2.287 


THE  TINES  IN  SEC  TO  BE  USED  FOR  INPUT 


TO  A  POSITION  ESTIMATOR 


TIME! 1 >= 
TIME<2>= 
T IME  <3> • 
T I ME  <  4  >  - 
T I ME  <  5 )  = 
T I  ME  <  A  >  - 


-5. 360  SEC 
-5.040  SEC 
-3.723  SEC 
-2.007  SEC 
.000  SEC 
2.287  SEC 


1  HI.  flET  TO  A  BP  USED  AS  INPUT  TO  POSITION  ESTIMATOR 


TEMP=  5.6  BEGS  C. 
WIND  DIRECT I  ON-  6380  MILS 
WIND  SPEEU=  a  KNOTS 


Figure  5 


Results  from  PASS  shot  S1-D3A1-MB-2 
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t»..j  *4i:  shot  rtm:  12:15  hrs  .mo  sec 

2  NIKI  ARRAY  B 

n  i  i,s :  3 / ,  ;*:j * 29 * 3u r  3  i  . 32  iv:.*o  val uu  s/mi c: 

TUI  -  ini;  (  INC  UMES  I'  01-  EACH  MiCRUPHUNE  TONAL 


f  i C  ( 

.'7 )  * 

12: IS  MRS 

35 

Civ 

SEC 

MIC  <  30) : 

12:1 

5  HRS 

y.j 

22  7 

SEC 

rijc< 

28 ) : 

12515  MRS 

3  5 

81  9 

si:  c 

n  i  c:  <  3 1  > : 

12:1 

5  HRS 

36 

SO  7 

sex 

M  1 C  ( 

29 ) : 

12515  HRS 

34 

331 

SEC 

N  J  C  < 32 ; : 

12 :  j 

5  HRS 

3i.J 

297 

see 

WINDOW 

STARTING  TIME 

3  (MSEC 

>  RELATIVE 

TO  A DOVE 

MIC  1 

MIC  2 

MIC  3 

MIC  4 

MIC  5 

MIC  6 

623 

733 

730 

746 

779 

701 

S  J  UNAL 

ENERGY  WITHIN 

EACH 

WINDOW 

nil:  i 

MIC  2 

MIC  3 

MIC  4 

MIC  5 

MIC  6 

2/5637 

.  35194 

377331 . 

540671 . 

330044 . 

362904 . 

1 1 U;  CORRELAT  TON*  1 1  ML  DIFFERE  NCFSUISE  C  )  r  AND  CORRECTED  TIME  DIE  I'ERCNCL'J SEC  » 

in: run  k  each  iatr  or  microphone  signals 


GAMMA (  1  *  j ) 

.  455 

K  t  AU  <  J  y  2  >  ■ 

14.0 

C I AU C  l  *  2  >  - 

-.014 

Mi'h  111  A  (  l*3>  • 

.  529 

R  r AU •  J  y 3 ;  • 

■  *490.0 

c:rou<  1 ,  5  >  = 

.  4  9f> 

gamma  <1*1  > 

♦  663 

l;  T  All  (  1  *  4  > 

- 1 399  *  0 

r  TY-iU  (  1  .  4  >~- 

•  1 .  ??■;• 

•  1 A  Mil  A  »'  .1  r  S ) 

.  420 

KTALK  1  f  5>:- 

•2708 . 0 

C  I  AU  (  1  r  5  )  - 

•2. 70 7 

GAMMA  •:  L  ?  6  1  - 

.  •’  1  0 

\i  T  AU  (1*6'  : 

4110.0 

C T  AU ( 1  .6 ) - 

-4.418 

GAMmA  <  2  *  5  .*  •- 

.  44  i 

KT  AU  (  2  r  3  >  ' 

4154  » 0 

CTALK2.3) 

-  .  484 

GAMMA  (  2  r  4  >  - 

.4  72 

KTALK  4  > 

1 384.0 

C  T  A  U  (  2  *  4  ' 

-  1 . 395 

uAilMA  (2y'jJ 

*  1  1  3 

R  TALI  (  2  v  S  '  - 

••  2692 . 0 

C T  AU ( 2  *  3  > 

-2.692 

<  j  AMMA  C  2  * 6  ) 

.  4 :  j 

K I AU *  2 »  6  > 

-4404.0 

r TALK  2* 6  i  = 

-  4 . 403 

GAMMA  <  3  *  4  .< 

.581 

R T AU v  ( t  4 : 

-901 .0 

C T  AU  <  3  *  4  » • 

.901 

GAMMA <  -• 

.  3  64 

I-’  f  AU  <  3  *  5  ) 

2206 . 0 

C  T  AU ( 3  t  5  > 

2 . 20° 

*  ■  AMMA  ( .i  *  6  *  : 

.  58  \ 

UTAH (3.6) 

3920.0 

CT  AU ( 3  y  6  > - 

•  3.9  3C 

GAMMA  (  '1  *  5  > 

.  4C2 

RTAIG  4fS>- 

J  308.0 

CTAU(  AfVj) 

: .  5on 

>  .AMMA  ( 4  f  6  > 

.626 

R  I A 1 1  (  4  y  6  ) 

-3018.0 

CTAM< 4  *6  ; 

3 . 0  j  9 

GAMMA  ( 5  t  6  ) 

.314 

IV  I  AIK  5 . 6  >  - 

1710.0 

C  T  AU (5*6) " 

- 1 .  /  n 

TIE  TIMES 

ru  si:  c 

TO  BE  USED  FOR  INPUT 

ru  A  POSITION  ESTIMATOR 

r  1.  Ml  (  1  ) 

-1.399  SEC 

1  I  ME  {  2  )  ~ 

-1.385  SEC 

r  .1  mi:  <  >)  =- 

*90]  SEC 

i  :r  mi  -  <  1  >  - 

.000  SEC 

TIME (5) 

\ . 30b  SEC 

rihi  <6) 

3.0J9  SEC 

Mil  Mi:  )  TU  A  DU  US),  n  A:  jNfUT  10  POSIT  If  IN  ESTIMATOR 

TEMP-  7.9  nr  OS  C. 

U  INI '  DIR  E  C  7  1  UN  -  AO  SO  MILS 
WIND  SI  'I  f  1 1  ••  8  KNOTS 


Figure  6 


Results  from  PASS  shot  S2-D341-MB-1 
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bay  31  i:  suor  iirtr  u:i 5  hrs  i  . o?a  sec: 

■JOUkU.  r  rtJKIL  aura/  k 

M  1 CS  :  27 *  2U  *  12V  V  .50  *  5  i  ,  3?  1*20  VALUES  /M  i  C 

t hi:  starting  times  roc  each  m  i crophune  signal 


ri  l  C  <  27  )  i 

0  « 1 8  HRS 

*4  .  V71  SEC 

MIC (30) : 

8: 18 

HRS 

'So .  379  SLC 

filC  (20 > : 

o:  ri  HKH 

34.971  SEC 

m  r  c;  (  3  i  )  : 

o :  l r. 

HRS 

37.689  ST C 

a  re  <27 ) : 

Otis  HRS 

38.483  SEC 

MIC < 32  > : 

8 :  i  s 

HRS 

39.481  SEC 

ri ir:  wjnuuw  s  tart t ng  t i meg r msec »  kela  i  roc  to  above 


me  1  MIC  2  MIC  3  MIC  4  MIC  3 

742  726  677  666  686 


THE  SIGNAL  ENERGY  WITH JR  EACH  WINDOW 

MIC  L  MIC  2  MIC  3  MIC  4  MIC  5  MIC  6 

143314V,  244964.  0322272.  21142428.  32418847.  34393879# 


r Hi:  cuRRELivnoNf  t j m  hif  iiirences ( msec .» *  ani»  corrected  iime  r« x r  Terences <oeo 

BETWEEN  EACH  PAIR  OF  MICROPHONE  SIGNALS 


GAMMA < l *?> 

♦  4?‘  i 

RTAU< 1 >  3 > - 

5*0 

C.T  AIM  1  *  2  > 

.014 

GAMMA (If 1 ) 

.  703 

KTAUdf  3) 

-422.0 

C i AU ( )  *  3  ) 

-.436 

GAMMA (1*4) 

-  .669 

R  f  AU  (1*4)  -■ 

•  1324.0 

C  l  AIJ  <  J  *  4  ) 

1.340 

GAMMA ( 1 v  8 i 

•  68  5 

R  1  AU  (.1*5) 

-2652 ♦ 0 
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«  768 

R  r AU ( J *6> - 

-4  329.0 

C I AU (1*6) 

-4.319 

GAMMA  <  .’*  >> 
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RT  AU  (7*3'  " 
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l.  I  AU  (  2  *  3  * 

.451 

GAMMA < 2  *  4 ) " 

.  703 

U  ) AU ( 2  *  4  •  - 

-  .1  359 . 0 
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-  L  .  384 

GAMMA <  2  *  5 ) 

.  rso 
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-  2.656 

G/ irtfl  A  *.  2  *  6  ) 
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k TAG <2,6  » 
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C 1 AU 1  8  *  6 > 
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R TAG  <  A *4  >  - 
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C  I  AU  ■:  3*4) 
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i:  1  AU  \  >  *  5  > 
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GAMMA »  3  •  A  ) : 

.'734 

R  r  AU ( 3 • A  > - 

-3  i  /  7 . 0 

CTAIU  A*. .  j 

-  3 . 882 

GAMMA ( 4  ?  8  > - 

.  798 

\\  TAG  i  4  ?  5  >  - 

-  1298.0 

C  TAIU4 ,5  > 

-  i  .  303 
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,  753 
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TIME  t  4  >  ~  -2.  *,7?  SEC 
T  I  M  E  ( 5  >  -  ]  .  6  1 6  S  E  C 

T 1  ME  <  6  )  ---  .  000  SEC 


mi:  Mr. r  to  a  ri:  user  as  jnpi.ii  ro  position  fmjmmIor 
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Figure  7.  Results  from  PASS  shot  S2-D341-MB-2 
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II MEG  IN  SEC 

TO  Bfc  USED 
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TIME  <  1 ) s:  - 
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TIME (3)-  - 
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Figure  8.  Results  from  PASS  shot  S2-D341-MB-4 
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me.  <30>:  a: 

3  HRS 

34.843  SrC 
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.  Results  from  PASS  shot  S3-D341-MB-2 . 
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R  I  AU  CHS) 

391 7.0 

CTAU(2*8) 

< .  929 

unflMA ( 2 * 6  )  ■ 

.  5  L  7 

R  !  ALJ  (  7  *  6  >  = 

461 9.0 

C  i  Al  1  ( .  * 6 '»  - 

4 . 674 

Cm'. Mi  i  A  (  3  »  4  )  - 

,  sen 

RTAUC3.4) - 

1314.0 

r  i  ai  i  ( 3  *  m  ■ 

i  .2:  3 

GAMMA <  3 » 8 )  - 

.  SOS 

R  T  i-ill  (  3  *  V* ) 

73S3.0 

C  1  AIJ  (  a  »  S  • 

.  3 1 : 

f  mVIMA  (3*6)  = 

.  1 3  8 

Rf  A  U  (  3  *  6  N 

302  7 . 0 

C  I  AU  ( .■>  *  6  > 

3.0  36 

A  A  f  IMA  (  4  r  S  '  - 

.  G  A  2 

RTAI.K  4*8) 

1 004.0 

e  i  AU 

1 . 0  t  8 

GAMMA (4*6) 

7.03 

RT AU (4*6) - 

1/07.0 

r 1 AU <  4  •  A > 

1  .  •  i  3 

GAilhA  ( S  9  6  )  •- 

.881 

RI  AU  (  8  *  6 ) 

699.0 

l  I A  U  (  8  *  6  ) 

.  <'7-’  8 

TIMES  IN  see  TO 

HE  USLii 

FOR  INPUT  TO 

A  POSITION 

11  ST  1 MA  1  UR 

rinru)  4.781.,  see 
T  I  ME (2)^  2.VJ1  see 

r  [MI  (3)  1.323  sec 

T  1  ML  i  A  )  .000  St:  c 

TlME‘l5>-  1  .01  o  SI  C 

n  ml <&>==-  - 1 .713  see 


(Hi  her  ro  A  ftE  USED  A>  INPUT  10  POSH  JON  ESTIMATOR 

TEMP--  G.4  DLGS  C . 
Ml  Nil  DIRECT  I  ON  ,.100  MILS 
W  l  NO  SI  T:: LI'  V  KNOTS 


Figure  10,  Results  from  PASS  shot  S7-D341-MB-2 
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Table  3.  Miss  Distances  for  PASS  Data  Examples 


File  No. 

Fig. 

No. 

Miss  Distance 
(Meters) 

Distance  to 
Center  of  Array 
(Meters) 

Ratio  of  Miss 
Distance  to 
Total  Distance 

S1-D341-MB-1. 

4 

430.0 

12,678.5 

.034 

S1-D341-MB-2. 

5 

195.5 

12,678.5 

.015 

S2-D341-MB-1. 

6 

176.4 

11,784.8 

.015 

S2-D341-MB-2 

7 

192.8 

11,784.8 

.016 

S2-D341-MB-4 

8 

!  1015.2 

| 

11,784.8 

.086 

S3-D341-MB-2. 

9 

395.54 

11,471.9 

.034 

S7-D341-MB-2. 

10 

291.8 

16,900.7 

.017 
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V.  CONCLUSIONS 

The  present  edition  of  the  "Time  Difference  estimator  Program"  has  been 
presented  in  the  first  three  sections  of  this  report  and  the  results  from  its 
application  to  a  number  of  arbitrarily  selected  signals  from  the  PASS  experi¬ 
ment  was  given  in  the  fourth  section.  The  results  show  good  time  duration  esti¬ 
mates,  based  on  the  concept  of  miss  distance;  iiowever,  one  should  realize  that 
with  only  seven  samples  not  much  can  be  said  of  the  programs  statistical  per¬ 
formance  other  than  it  is  looking  very  promising. 

From  Table  2  and  Figures  4-10  it  is  seen  that  the  program  satisfactorily 
determined  relative  times  such  that  target  position  could  be  estimated  with 
miss  distances  around  two  hundred  meters  for  targets  of  about  twelve  kilometers 
in  range.  It  appears  that  as  one  would  intuitively  expect,  correlation  coef¬ 
ficients  much  lower  than  0.5  for  the  time  differences,  Figure  4  and  Figure  8, 
result  in  the  considerable  sized  miss  distances  given  in  Table  2. 

The  application  of  the  progiam  to  existing  data  has  pointed  out  need  for 
further  research  into  several  areas  in  order  to  improve  the  procedure. 

(1)  In  order  to  eliminate  inaccurate  data  a  threshold  value  for  the 
correlation  coefficient  must  be  determined  for  which  the  associated 
time  difference  estimates  with  less  than  that  threshold  are  claimed 
unreliable  and  discarded. 

(2)  A  procedure  should  be  developed  to  assign  realistic  weights  for  use 

in  the  least  squares  procedure  perhaps  as  a  function  of  the  normalized 
correlation  coefficient. 

(3)  The  overall  effectiveness  of  the  procedure  should  be  established  by 
providing  the  variance  or  a  bound  on  the  variances  of  the  time 


Jill  erence 
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estimates.  This  will  require  a  theoretical  development 
■‘>U,U'  wiLh  an  examination  of  a  large  amount  of  data  for  experi¬ 
mental  verification. 

Ihe  program  should  he  altered  to  present  a  workable  procedure  for 
handling  the  multiple  target  problem. 


I 
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APPENDIX  A 


Time  Differences  Estimator  Program  Listing 


rim  [i ii  r i 'ur.rji a  .  i  . i  i or, i uk  i h/I/KApi 

■  iLNI-kae  I'l  ,okipiiun:  1 1 1 r»  *  utji iKr.n  un  i  .-.f.r  i -f •  i  a  givin  nomi-i  r 
01-  M  ILRul  HUNK  S  I  GOALS  l  ACU  EONTAIUliJt  A  bin'll  J  1_  l  r  ■  •  *  - 1  .  i  : 

or  l*a i i'l  points  i  kun  a  i  i>lp  uni  a  <  10-1-1..  i  anu  hi  ni  n.o,  1  s a  r 1 01 

Hll  I  I:  hi  NILS  M:  TUI.  IN  SIGNAL  PURGES  I  HA  1  API  TAP  U  I  i  I  i  i  i «  Till 
RLl.OKl*  LENGTH. 


AUTHOR:  I.  (INN  il .  C  .  LIU'I:  MAN 

elec .  anp  coopuilr  eng,  iiij  i. 
ni:u  hi.  a  1  co  0  r a  r  1:  uk  1  vi  s  i  r  r 

l  AO  CRUEL G*  Nl.  U  hi.  XI  CO  UU003 

r  1: 1. 1;  phone :  sos-64*  1321 


GLOSSARY  OF  VAR  t  AtL  l“S  X 

NMJC-  NUMbE'R  01  0 I  CROPHONIiS 

NL~  NUMBER  OF  HAT  A  POINTS  I  UR  II  At.  I-!  O  J  l.TUPHUNE 
NSW  NUMBER  OF  I  '0  IN  TO  IN  SI.  1 1U  NO  UJNl'UU 

NPFW  -  NUMBER  01-  POINT''.  IN  rill-  I  uURJ  I.'K  UlNJiQU 

Si  ANA  I  AMPLITUiil.  FOP  0  L'C  J  A  I  llol  J  kl.l  AT  ILT  TO  IT  (  D 

iT(.n-  s  i  Apr  j  no  time  rop  n  i  irui-iione  j  in  iisi.  l  relative  n  pi  as. 

LRSMJ)--  ROUGH  SIARMNG  POINT  TOP  0  I.  CROFMH  'Nlf  .1  PI...L  A  T  J  Cl.  T '  •  II  <  U 

NR  (  .1  )  -  ".TART  T  l  ML  I  UP  I'UlIKIIr'R  WINf'UW  O.  n.lC  I  KELi.MV!  In  Til: 
E<J#  -  I:  Nil:  RO  Y  WITHIN  GORKI  'L.ATION  UINtsi.J  i  ()|>  Si  ON  At  ; 

TAIKIrJ)  RAW  T  T  HI  LU  I  I  TRI.rJCi:  Plr  TUIIXN  W  I  N; JE  H  S  I  GNAI  ’  J  ANG  I 

KTAUitpJ)-  ROUGH  1 1  Ml..  I*  Tl' F  L  FT.  NQ  M  l  l-'L  I.  N  r,  i .  keJ  ') IUNL  S  )  A  N I ' 

01.  I  t.Ri  ilNl'.ti  0  r  CORKI.I.  A  I  ION 

CTolJtlrJ)  lOPPHUHi  Tint  MTITRiNPI  I<l*  I  UEI' f.  n  U.ROF'MONi:  I  AN!'  J  FT 
LEAST  s  UUAIvF.il  ERROR  I  I  T  (  KL.I.  AT  10  L  * 
liAh  (!*.)>'  NORMAL 12EU  DlkRF.L  A  \  I  ON  I'JI  I  EFHTFN  hll  .1  i  SIGNAL  s 

Rlf  I  [AH  111  TV  UL  I  OUT  ASSIGNED  TO  KiAiU'ni)  TOR  I  HI. 

LEAS  I  SUUARi:  I  'ROC II  UlJRE 

A  T  <  J  >  llOl  Min' ERL:  NCI  Ul  SIGNAL  J  RLlATIVt  TO  hi'  i  SiuNAl 

I  I  ML  <  1  '»  RIL  ATILT:  TIM.  OF  THE  1  TU  GI'-NA..  I  OR  PuS  .1  T  1 1 1  J  r'Tihc'Gi- 

A<J.i;  AUGMENT  II'  r.OI  I  r  JCJtNT  h  A  I  R  l  X  I  UK  I  i.ASr  SUUARI  1  M 

»L  (  I  #  ■  (.  ON  f  KM  UAL  I  AIM  L  S  f  OR  SUORUU  T  1  Nil  » ■  IK 

LM  )  OAK  1 AJM. I  0  Vv  I  >  ANI*  V(2>  I  OR  INPUT  Tu  G  JR  GUBkGU  T  1  O’; 

ini  AS(I»-  OLM  KALI  hi  ASUIVL  Or  CORRELATION  FOR  h!C  1 

ill  CilbU  I-  MAKilNG  AN1I  MOPPING  FREUUENCY  V.’r  UF.G  F  OR  60  112 
COOP  ANH  LOU  PASS  V  II  ILK 

*.  (I>  J  TH  UAL  UL  01  FIRST  SIGNAL  FUR  CROSS  COR  RE!  LA  I J  ON-  INI-U) 

I  Til  UAL Uf  01  I  OUR  !  I: R  TRANSFORM  UN  OUTPUI 
Y  <  J  )  1  HI  UAl.Ul:  Ul  SI-  CONI1  SIGNAL  L  UP  CROSS  CORRELATION  INPUT 

I  TH  VALUE  UF  FUURJER  TRAN SI  UPh  ON  OUTPUT 

UOK’K  «.  I  >  -  UGLTv  .  MIL  SPLf:  .1  L  T  LU  E<  Y  I'UURG  SUlMcOUTtNI. 

xm  i  )  j  iH  value  oi  rsTihon:  of  noisi  i.  nergy  si-ll  ti  on 

Y  Y  (  T  )  1  TH  VAL.UL  01  I  ST  1  On  1 L  Ul  SIGNAL  PLUS  NO  1ST  SlLlTRUh 

U.-.l  Ul  or  PRFIMI  UK  AT  Tl-ll  I  !  II  I  Rl  Uiu  NCY 
TIiAYiJi  Li  A  Y  01  THE  START  01  T  HI  I  IH  OICRUPUONE  S1C-NAL 
1  HOUR  <  f  >  HOUR  01  HIP  START  OF  TNI.  3  111  hUROPHONI  SlUNAI 

]OlN<r>  ri  I  NT  IT  E  01  TNI  SIAKI  01  THE  ITU  r :  ]  I  ROPI  TONE  SIGNAL 

SI. IS.  1)  SE  CONI'  01  THE  -ImKI  01  THi.  i  OI  1 1 1 1  R(M-T  TONI  S  I GNAL 

Oi:C(jM  I  ASS  OICROPIICINI  NHNEI.R  ASSIH  lAiri'  WITH  hji  ROPHONE  l 
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.vjJi  INI  I  tAL  12  INC, 


61 

A  1 

c: 

U I  Ml.  NO  1  (JN  1 H A  Y  (  A  )  »  1 1 IOUR  <  A  )  »  J  M  t  N  <  A  >  » 

SC  0(6/  »  M  1  C  (  A  )  p  T  I  ME 

63 

H  1  ML NS  I ON  1  Al  1  f  A  »oJ  »  J  1  (  A  )  *  NR  <  A  >  r  F  ( A  > 

r  \<  1  AU  ( A  f  6  >  *  C  I  All  (  A  »  6 

6  4 

Ii  1  ML NS  1  ON  Sh  I C  *  A  *  1970)  r  f  RG  F  <  6  >  *  G  <  A  » 

A  )  r  A  T  <  A  )  ,  A  <  5  f  A  >  r  Jl.  ( 

Ali 

£i  1  Mi:  N; J  J  ON  CMEfiS  i  A )  .  hCUMK  (11* 

A  6 

CUMFLCX  X<  1024  ;  t  Y(  1  024  >  f  WORK  1  074  *  t 

XX (  1074  > .  fY<  107  4  >  »2 

A  7 

NM.U  A 

60 

NL  •  1 970 

69 

NSW- 347 

70 

Ni  l  U  5  .12 

71 

c 

72 

c 

J&R 

INGING  IN  DATA  FROM  MASTER  HU 

73 

c 

74 

CALL  RHAI  A<SnlL»rt.[Cr  lHAt'*  IHUlJK.  LrUN 

f SEC f TEMP, MILS* KNOT 

75 

DO  50  I-I-NMIC 

76 

I  SEC  =  SL:  C  >  I  >4  1000 

77 

NU  >  ~JIIOUR<  I )  *60460*1000 F  CM]  N*60t 10004 ISEC 

78 

50 

CONTINUE 

7  V 

c 

80 

i: 

r 

COMPUTATION  OF  ROUGH  ARRIVAL  UrttS 

o  1 

82 

DU  100  J-'lrNMlC. 

83 

1 E  =  0 . 0 

84 

t MAX -0.0 

88 

LMAX -0 

86 

DO  90  l\  1  *  NSW 

3  7 

i;  “Sri  n :  (  J  *  R  >  4>.f i  T  C  i  J » K  > 

33 

TE  - TL  fl. 

89 

90 

CONTINUE 

90 

EMAX " TE 

91 

LMAX-0 

9  8 

NLr’-NL  -NSW 

93 

DO  95  l\!:s  1  *  NLP 

94 

T  E  •*  TE  +  Sri  t  C  <  J  »  NSW  { 1*. )  #5M1C  *.Jr  NSW  *  R  >  - 

sri  if.  <  j< i,  •  »  ;,m  tc t  j*i- 

98 

r r  c  n_ .  lt  .  ema x  ,•  go  to  95 

96 

F MAX "T E 

97 

LiiAX-K 

98 

95 

CONTINUE 

99 

l.RST  ( J  >  -L  rlAX 

LOO 

100 

CONTINUE 

101 

c: 

102 

c 

COMPUTATION  Of  UINUOU  STARTING  POINTS 

.1.03 

c 

10  4 

DO  110  1-irNMIC 

105 

I NSW--NSU/4 

1  06 

MRS  I  •••  L.K’S  T  <  I  >  I  NSW 

107 

IFCMRST.GT .0)  GO  TO  105 

108 

NRO  >“0 

1 09 

GO  TO  110 

1 10 

J  1 1 

105 

CONTINUE 

NP:--NL-NPFU 

J  1.2 

IF  (MRvr.L  r  *NP>  GO  TO  1.06 

1  I  .? 

NR  <  T  >  tNI  * 

J  14 
i  1 5 

I  T  (  1  )  l  T  l  '•  4  NP 

GO  TO  tiO 

1  LA 

•  06 

ii|<(  .1)1  1  I  <  l  >  -CN:,U 

1  17 

T  T (  1  >  JT< 1 >+NR( 1 > 

1 18 

1  I  9 

u 

1  1  0 

CUNT  INUI 

f  O  A  i  i  (  f>  f  A  > 
>  »  V  (  2  ) 

<  i  o:m  > 


Ik  lil  TING  Till.  Ul.UiUU  i  I  .MV  |  IMG  Ufil-S 


li'jJ 
1 36 : 
137J 
138: 

139: 

140: 

1 4i :  c 
I42:c 
J  4  3 :  c 
144: 

1 45:  ; 

146: 

147:  ; 

148: 
149: 
100 :  . 

iv,  i  :c 

v. j2  :c 

r:3  :l 
!'•  4; 
100: 

1 06 : 
to?: 

1 08 : 
109: 
160: 
i6i : 
I62: 
163: 
164: 
160 : 
166: 
167: 
16Q : 
16?: 
170: 

1 7 1 : 

1  72: 

1 73 :  c 
174  :c 
i/ojc: 
1/6: 
177: 

1 78 : 
179: 


<41- 1 M  <6, 30 11- 

'Oil  ruKIt,-. !  (  /  /  '  Ilil  WINDOW  S  1  AI-:  I  I  NO  I  I  fil  5(M!.L 1.  -•  KIN  At  ILL  TIJ  AD'JVI 
UF,  III  ( t.  •  .'0  L  3  ) 

t.<  I  DKM.'i  I  (/ '  MIC.  I  MFC  7  MIC  3  M  l  L  4  dll 

.  II.  6'  '  > 

UK- 1 1 11  <  s,  :;•./}  o>nr 
.-010  em.-ma r ( ix .<■><  1 10)  > 

comfijiai kin  ui  signal  i.nfkgy  within  fourier  window 

Ml  1  15  J  1  .NMIC 

SUM  0.0 

RSI  AR  NR(J>+1 

KOTOR  -NR C J >  +NF'FW 

Ml  11.1  l\  - 1\  5  f  A  K .  I\LT0 I 

f=SlllC(  J,K)*SMIC:<  J.K> 

SUn ; SUM  TP 
1  1 2  CONFJNUL' 

115  E ( J )  SUM 


PRINTING  lIC  SIGNAL  ENERGY  UJIHJN  LAUI  WINDOW 
UR [ IE ( A  r LOO ) 1 

joo?  format  <//'  nn  signal  lnlkoy  wiihin  i  acfi  window'/i 
ur-:i  rt<-Sr20i4) 

101-1  F  ORMAT  C  '  MIC  1  MIC  .7  Mill  3  MIL 

MIC  5  MIC  6/> 

WRITE  IA»200!3>E 
..’003  F  OkMA  I  C  6  (  FT  3 . 0  )  > 

r-O'-n  r  I  OWING  OF  SIGNAL:-  I  OR  DETERMINING  1FIL  CROSS  C0RRF1LAI  ION 

DC)  .’00  l  lr'j 
M  =1  +  1 

M)  ISO  J  M  .NUJC 
RA-NRI I ) 

DO  120  I  =1.517 
KAL=RA+L 
SN-=SMICn  .RAL  ) 

X(L>  --CMF-L  X  ( SO.  0 . 0  > 
y(i..)--m:mpi.x  ;o.  o-o.o) 

120  CONTINUE 

DC)  130  l  =513.102-1 
X  <  L.  )  -CMF'LX  CO. 0.0.0) 

KK=NR<  J) 

KA'.!-=KBIL-r.l2 
SKAL=SMIC< J.KAl  P> 

Y  C  L  >  “ChF'LX  ( SKA!  »0.0> 

130  CONTINUE 

CALI..  I  UURG  <  X  .  1 071  >  -  t .  WORK  ) 

CALL  F  OUR'G  ■.  Y  .  1024  .  -  1  .  WORK  ) 

DC  TERM (NAT  ION  OF  STATING  I  'UJNT  FOR'  NO  1 : -  F  WINDOW 

KXX-NR ( I )  -5  1 3 
KYY-NR'C  J  1  -51  3 
DO  125  L  1.512 
I  l  1.  +  517 
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i no :  xx (i.  i  >  . c . hi  i  >.(o.o 7 v» . o > 

i ;  t  [ :  rr<Li)  Ctin 

USj:  ll  Cf\>  X  *1.  !  .  J  ;Ml  IU  11A 

1U3 :  mxx  =  i*  >; >.  i  i. 

I'JAt  3  XX  Sri 1  f J  C  I  *  MX/  ; 

t  u *:*  i  x X  <  L  >  r nr-'j . aO-xx»  o  .  o  > 

Uf6J  i  10  CONTINUE 

UU:  If  IKY  V  *L Y  •  1  >t»0  10  117 

liK<:  myy~k  yy+l 

t : )  9 :  s  Y  Y  --SM  X  C  (  J  f  H  r  Y  .) 

1  ••'o :  y  r  <  l  >  chpl x  (  s  y  y  r  o .  o  > 

1 7i :  go  ro  1 25 

1-72:  Uo  XX  <  L  )  : CMPLX  <  0  *  0  *  0  *  0  > 

vn:  go  ro  uo 

I'Mi  H?  Y  Y  <  L  )  -’Cfir-'U  X(0. 0*0.0) 

175 S  125  CONTINUE 
1 7  6  :  c 

197  :c  DEIERMINAf  ION  0/  C5TIMA  Hr.  OF  NOISE  SPECTRUM 
190  :c 

3  99J  CAL  L  FQURG (XX 1 1024 »  - 1 » WORK ) 

200  :  CALI.  FQURG  (  Y  Y  ,X  024  *  -  1  ■■  WORK  ) 

201 5  DO  135  V.--1  f  1024 

202 :  ZZ< L  > -CON JG < XX  < L ) ) 

20 3 .*  XX <L  >  fXX(L)*ZZ(L  ) 

20  A  :  ZZ(  L  )  -■  C  ON  J  G  (  Y  Y  ( L  )  ) 

205J  YY(D  ~YY<L>*2Z(L> 

20 o  :  Y  Y  <  L  >  *  ( X X  ( I.  )  +2Z  <  I... )  ) /2 . 

2o.':c 

208  tC  DETERMINATION  C.i\-  ESTIMATE  01-  SIGNAL  PL  US  NOISE  SPECTRUM 
209 :  i 

21 o:  ZZ( L>  Y<L> 

21 1?  Y  <  L  )  'CCjN  JG  <  Y  (  L  )  ) 

212:  X  X  <  L  >  ~  Y  (L  )*ZZ(L  ) 

21  i  :  27  (I.  >  -  CON  JG  (X(IJ) 

714:  Z2<L>  ?T<L)¥XU.  > 

2  I’jI  XX  <  L  >  -  <  XX I L  >  +77  U ... )  Y  / 2  , 

:ma:c 

7 1 ; :  c  compu  r a  r  ion  of  pro  f  ltek  ual u>:  s 
218  :c 

219 :  PKA-*REAL(XX(L)  ) 

220  :  FRB-REAl  O  Y(L)) 

22 1:  PRTi-.OOl 

2  2  2  :  IF  ( PR  A .  GT  *  PKTi )  GO  TO  7 1 6 

22 T:  FRFF-l.O 

224.*  GO  TO  721 

225 J  71 A  CONTINUE 

226  :  PROPRR'PRA 

22?:  ;r f  <  prc  •  gt . a  » >  go  to  7 1  v 

228 :  FRCP *  1 . O-PRC 

279 t  GO  TO  721 

230  I  719  FREE  =0,0 

231 J  72 A  CONTINUE 

2.52 :  ZZ( L  )  =CMI  LX  < PREP  *  FREE  ) 

233.*  Y(L>  -X(L  >*Y'L> 

234:  Y(LUY<L>  /.l 024, 

735 :  C 

236  50  APPL1CA1  ION  Of  mni  ll  K 
23  7 :  C 
230 : 

23  9 :  c 


Y  <  L  >  -  Y  <  L. )  *  £  2  <  I  > 
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APIL  t  CAT  ION  01  COMB  I  II  N  R  FOR  60  CYCLE  AND  HARMONICS 


135  CONIINUL 

Bi  i  I  A  MCOMB/1  r-  *  10 1  60  r  Aj  »  121  r 126f  J83f  1  HU , 244  1 2  4Y/ 

N  *  •■•(iCUMb  <  1  ) 

II  (NX.LE.O)  GO  TO  131 

Y  <  1  )  -CMPLX ( 0 . 0 » 0  . 0 > 

131  DO  133  N:=2f  10*2 
NX-MCOMB  <  N  > 

IF (NX.LE.O*  GO  TO  138 
NY  "MCOMB<iM  f  1 ) 

DO  132  L I.  =NX  »NY 

Y  <  LL )  CMPLX  ( O « 0  v  0  •  0  ) 

MX  2INI:'I::  W  K?  LL 

Y  <  MX  >  -CrtPL  X  <  0 . 0 f 0 . 0  > 

132  continue: 

133  CONTINUE 
138  CONTINUE 

FINDING  THE  RAW  TIME  DIFFERENCES  BETWEEN  THE  WINDOWED  CIGNA IS 

CALL  FOURGCY*  1024 f+1  rWORK) 

CMAX=0. 

TMAX  0 . 

DO  MO  L‘- 1  v  1024 
LA  REAL  <  Y  <  L  >  ) 

CA -ABG ( CA ) 

I F  <  C A « LE • CMAX )  GO  TO  140 
CM AX  CA 
TMAX-L 
140  CONTINUE 

TAU <  1 » J )  TMAX  NF’FW- 1 

DETERMINATION  OF  THE  NORMALIZED  CORRELATION  COEFFICIENT 

G  < 1 1 J ) -CMAX/ ( (E  <  J  >  *E ( I > ) **0 . 5 ) 

GAM iI»J)«G(IvJ) 

ISO  CONTINUE 
200  CONTINUE 

DETERMINATION  OF  ROUGH  TIME  DIFFERENCES  BETWEEN  SIGN  ALL.  I  AND  L 

WRITE (6? 20 16) 

WRITE <6 *2017) 

2017  FORMAT ('  BETWEEN  EACH  PAIR  OF  MICROPHONE  SIGNALS'/) 

DO  250  I«1»S 
KX- 1+1 

DO  240  L-KXrNMIC 

RTAU ( I  * L  ) -FLOAT  ( 1 1  (  I )  -  IT  ( L  )  >  +  T  All  (  T  * L  ) 

2016  FORMAT (// '  THE  CORRELATION*  TIME  DIFFERENCES (MSEC ) *  AND  CORRECTED 
.  TIME  D IFFERLNCES ( SEC  >  '  ) 

240  CONTINUE 
250  CONTINUE 

SELECTION  OF  MIC  TO  BE  USED  AS  REFERENCE 

CMEAS  (  1  )  G  (  I  *  2 )  4G  (  i  r  3  >  I G  t  1  *  4  '  f  G  <  1  r  5  >  f  (3  i  1  *  6  > 

CMLA3 <  2 )  G ( J  *  2  > I G ( 2  *  S) f  G ( 2 *  4  » LG ( 2 r 5 ) f G ( 2  * 6  > 

CM  LA  5  ( 3  >  -  •  G  <  1  ?  3  >  +G  <  2  »  3 )  K»  (  3  »  4  >  4  0  ( 3  •  5  >  +G  ( 3  •  .s  » 

CMEAS  <  A  )  *  G  ( 1  f  4  )  K.  ( 2  f  4  )  L C.  ( 3  •  4  )  LG  <  4  f  5  >  Hi  <  4  *  6  > 
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5'!.*:  cm  u<  i  .g>  hit  j»5>  * c; ( 3 , 5 >  haa*'.,)  h;(‘jfA» 

3-j  1  :  I  .ill  .V..  (  A  )  l>  (  1  •  /  )  K.  ( J  *  A  >  f  G  (  5 » A )  f  (j  (  \  »  6  )  H I  <  3 »  *  ) 

40 u.  ■.»  «  *.) 

4:  r«* i  .",\j  r  i  »6 

.404:  n  (i;mi_«v»<  j>  .l  r  .cc  >  go  to  233 

u.  i.rti'.MSd) 

io.*,:  k  - 1 

.40/:  KrtA/  k 

.so o:  255  lonunue 

.5ov :  r 

jio:c  i.stadl  tsumiNT  or  weights  for  least  square  time  i  it 
:>  1 1 :  c 

it.:':  do  2/3  i  i.'j 

iii:  Ku-iti 

.514:  DO  2/0  .  J  -  I\IJ  ?  A 

31 G:  IF <G(If J) .L! .0.5)  GO  TO  2AA 

3  I  A :  2  AS  G  <  1 » J  >  -GAM  (  l  y  J  )  *  1 0 . 3 

31/:  GO  TO  2/0 

318 :  266  G(  I,.l)-‘GAfi(  I  » J)*GAi1<  1»  J) 

319:  2/0  CONTINUE 

32 0:  2/5  CONTINUE 

32  l  SC 

3225C  DETERMINATION  OF  II.  AS  I  SQUARES  TIME  I"  IT  WITH  MIC(l)  AS  REFERENCE 
323  :c 

3241  HO  290  I -Ip  4 

220 S  KX-I+1 

32A:  DO  280  J-lvX  y5 

32/:  (Silt  J>  -G<Ifl  v  .14  1  ) 

328S  A  <  J  y  I  ) -A ( I r  J ) 

329:  280  CONTINUE 

330:  290  CONTINUE 

3  3 1  :  A  ( 1  y  J  ) :::  G  <  1  y 2  )  TG  <2p3)46<2f4MG(2*  5 >  +G  <  2  y  6 ) 

332  :  A  <  2  f  2  >  --G  <  1  f  3 )  +G  <2f3)4G<3f4>  EG  (  3  *  5 >  f  G  ( 3*  6 ) 

333:  A ( 3  r 3 ) :~G  < 1  *  4 ) f  G  <  2  r 4 ) T  G ( 3  t 4  >  + G ( 4  v  5 ) +  G ( 4  *  A ) 

334  :  A(  4  f 4  ) ;;-G  ( 1  y  5  )  +  G <2 y 5  >  Ki  <  3f5  )  f G  (  4  y  5 )  TG  (  3  p  6 ) 

335 J  A  <  5  f  5 ) -G  < 1 f 6 >  +G <  2  y  A)f  G ( 3 «  A ) f  G  (  4  y  A  ) 4 G ( 5 p A  > 

336  :  AA;=~G  (  1  y  2 )  *  TAU  ( 1  f  2 )  +G  ( 2  y  3  )  *TAU  ( 2  f  3  H  G  ( 2  r  4  )  *  T  AU  <  2  y  4  > 

3 37  J  A  ( 1  y  6  ) ~AA  FG  <  2  f  5 )  *TAU  <  2  •  5  >  EG  (  2  y  A  )  *TAU  <  2  f  A  > 

338  t  E<B  -  -G  <  1  y  3  >  *  T  All  <  1  f  3  >  ~G <  2  f  3 )  *T  AU  <  2  •  3  >  +G  <  3  f  rJ >  *TAU  <  3  y  5  > 

339  J  A <  2  r  6 ) «BB+G ( 3 f  4  >#T AU  <  3f 4 ) +G ( 3  p 6 >  * ! AU  <  3  f A ) 

340  J  CC”-G<  1  f  4  )  *TAU  1  t  f  4  )  C»  <  2  f  4  )  *  T  AU «  2  p  4  )  ~G  (3*4)  *  T  AU  (  3  y  4  ) 

34 1 J  A ( 3  y  A ) -CCFG ( 4  y  5 ) *T AU <  4  y  5 ) +G ( 4  y  A ) *  T  AU <  4  y  A ) 

342J  DD--G  <  1  y5)*TAU<  l  y5> ~G<2 y 5 >*TAU< 2p5 >~G C 3 r 5 > *T AU ( 3 p 5 > 

343  S  A ( 4  y  A ) “DD-G <  4  y  5  >  *T  AU  <  4  p 5 ) FG ( 5 y  A ) *T AU ( 5 y  A ) 

344 :  EE--G <1 1 6 ) *TAU < 1 1 6 ) -G ( 2y  6 ) *TA U (2r6> -G (3 1 6 M TAU (3r6> 

3  45  J  A  <  5  y  A  )  =EE-G  <  4  r  6  )  *  T  At)  <  4  y  A )  -  G  (  5  y  A  )  *TAU  (  5  r  A  ) 

34 A:  0(1) -A 

34/:  CALL  GJR<AyAy5y5yAyi400y  .JCy  V) 

348?  A T ( 1 ) ~ 0  « 

3  49:  DO  295  K- 2 y A 

350  J  AT  <  K )  :=A  ( K-  1  f  A ) 

351 J  295  CONTINUE 
352  :c: 

353 :c  ADJUSTMENT  OF  TIME  D IF  I  CRENCES 
354  :c 

335:  DO  310  T " 1 f  5 

336  S  KT-T  +  t 

35/:  DO  300  J-  R  r  r 

338 :  r au < i f  j  >  a  r  1 1 > -  a  t ( j ) 

359  :  c  r  AU  ( I  y  J )  - 1  AU  <  I  y  J  >  H  I..  OA  T  ( I  l  ( I )  -IT  <  J  >  ) 


Cl  All*  1  y  J)  -t:iAU(  1  r.D/1000 


c:  i  RJNriwu  cuKkELA r ion  cun  and  rough  and  currcc n:n  iim_  nin  pklnci.s 

c 

WKM  I E  <  A  y  200  l  )  I  »  .JyGAM<  I f  J ) r  3  f J r R I AU < 1 »  J) »  T  f  I »  C  I  AIJ(  (  >.J) 

200  L  FORMAT  (  AX  y  '  GAMMA  <  '  2  1  »  '  t  '  .7  J  r  '  )  'FA  RT  AU  ('  l  .1  r  '  r  '  1 1  *  ‘  )  '  f  /,  J  y  ./ 

.  X  y  '  C  TuU L 1  v  '  >  -  '  I-  A  .  3  ) 

300  C0NI1NUE 
3 J 0  CONTINUE 
C 

C  BETERMINAT  tON  OF  RELATIVE  TIMES  FOR  MIL  POSITION  ESTIMATOR 

c: 

no  A 20  I-lr5 

KX-I+1 

no  A  to  L-KX  *  Nil  1C 

CTAU  <  L  • !  >--CTAU<IfL.) 

610  CONTINUE 
620  CONTINUE 

HO  640  J  - 1 »  6 
IF < J«NE • KM AX ) GO  TO  630 
T  JME  <  J ) -0.00 
GO  TO  640 
630  CONTINUE 

T I ME <  J > -  C!AU< J.KMAX) 

640  CONTINUE 
C 

C  PRINTING  THE  RELATIVE  TIME'S  FOR  POSITION  ESTIMATOR 
C 

WRI  T E(6  y  2005) 

20015  FORMAT (//'  THE  TIMES  IN  SEC  TO  BE  USED  FOR  INPUT  TO  A  POSITION  ESI 
.  IMA  TOR  '/> 

HO  503  I 1  r  NMIC 

WR I TE  <  6 f 2003 ) T  f  TIME ( I ) 

2003  FORMAT  <  30  X  f  '  T I  ME  (  '  1 1  y  '  )  :=  '  F  7 . 3  f  '  SEC  '  ) 

503  CONTINUE 

UR I TE  <  6 y  2025 > 

2025  FORMAT (/ ) 

C 

r  OUTPUTTING  THE  MET  INFORMATION 
C 

URITE  <  6  f  2021 ) 

2021  FORMAT ('  THE  MET  TO  A  BE  USED  AS  INPUT  TO  POSITION  ESTIMATOR'/) 
URITE (6 f2022> TEMP 

2022  FORMAT <32Xf'  TEMP®  'F4.1f'  BEGS  C.  '  ) 

URITE <6 f2023) MILS 

2023  FORMAT (22X f '  WIND  DIRECTION®  '  14y  '  MILS') 

URITE  <  A y  2024 ) KNOTS 

2024  FORMAT ( 26X  r '  WIND  SPEED®  'I4y'  KNOTS') 

GO  TO  401 

400  URITE (Ay 2015) 

2015  FORMAT (30H  AN  OVERFLOW  HAS  BEEN  DETECTED  ) 

401  CONTINUE 


413 


Ei 
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APPENDIX  B 


Subroutine  RDATA  Program  Listing 


i  :c 

2!C  SUhROUMNE  NBA  l  A 

3  :c 

4!  SlIUKOUT  INIi  ETiAT  A  <  SMIC . M  I  >  I B  A  V  *  T  HOUR  •  IMIN.SEC.  nXrfUXLSrl'.NOI'j) 

5!  I1  l MENS I ON  MIC(6)  * 1BAY<6)  . 1  HOUR  ( 6  >  .  I  i-r  I N ( 4  >  .  5L0  ( 6  )  .SMIC<6.  1920) 

{>:  INTEGER  VALUESCA) .VAL 

/  !  KLAIK  12.  lll)MIMOT.MSORC*MB.MH*MMlN*SSEC*TEMI-',mLU*i  HOTS 

!3:  111  FURrtAI  <13. 13. IX. 13. IX. 12. IX.  12 . J  X >  F4 ,3»  3X.F  4  .  I  ,  1  X .  i  4,  IX  »  I  2  > 

9!  HH=hlL7 

10 :  WRIIE<4»115> 

1 .1  :  US  FORMAT  <///> 

121  MR I IE <  6 , 1 12 ) MV .  MH . Mrt I N . SSEC 

131  1.12  FORMAT  (4X* '  HAY  '13.'!  SHOT  TIME  'I2*'J'I2»'  HRS  '16,3.'  SEC') 

141  READ<12.221> 

151  22:1  FORMAT  (5X.50H  ) 

141  WRITE (4. 221 ) 

171  HO  11  I=1,MIKT0T 

10 1  REAR (12*333)  VALUES (I J , MIC ( I ) * IHAY< 1 ) . I HOUR < I )  *  IMIN<  1  >  *SCC.<  I ) 

191  I HOUR ( I >  =  I HOUR ( I ) -7 

20!  333  F0RMAT(3X»I4»5X. 13*14*13*13* F7 . 3) 

211  II  CONTINUE 

22!  UR  I  IE  (6  >225)  MIC  (1 )  *MXC(2)  »MIC<3)  >  MIC<  4  >  .  h  IC<  5 ) »  MIC.  <  4  >  *  VALUES  ( 1  > 

23!  225  FORMAT <4X*'  MIX'S!  '  12*  '  .  '  12*  '  *  '  12*  '  »  '  12*  '  *  '  12*  '  .  '  I2.3X*  14,  '  VALUES 

24!  ./MIC') 

25!  WRITE <  4*227 ) 

26 :  227  FORMAT!/'  THE  STARTING  TIMES  FOR  EACH  MICROPHONE  SIGNAL'/) 

27 J  HO  25  1-1.3 

20!  J  I ‘ 3 

291  WRITE  (4 .229)  MIC  (I ) ,  I  HOUR  ( I )  *  IMIN  (1 )  .SEC  ( I )  .MIC  <  J  >  » I  HOUR  <  ,J) , 

30!  . TMIN< J) ,SEC< J) 

31!  229  FORMAT <4X*'  MICl'12.')!  ' 12 . •  1  ' 1 2 * '  HRS  'FA. 3.'  SEC ' 5X . '  MIC('I2 

32!  ..')!  ' 12* ' 5 '12. '  HRS  'F6.3.'  SEC') 

33!  25  CONTINUE 

34!  no  22  I~J  .MIKTOT 

35!  VAL- VALUES ( I ) 

34!  REALM  12.444  >  ( SMICCI * J ) . J-l . VAL ) 

37!  22  CONTINUE 

38 !  444  FORMAT (16F5.0) 

39 !  RETURN 

40!  ENB 
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APPENDIX  D 

Derivation  of  Weighted  Least  Square  Solution 


The  performance  index  to  be  minimized  is  given  by 


IN  IN 

=  y  y  y.  . [r. .-(t.-t .) y 

i>j  j-i  1J  1J  1  1 


In  this  expression  T  represents  the  estimate  of  the  rougli  time  difference 
between  signals  from  microphone  i  and  j,  y _  is  the  weight  associated  with 
that  estimate,  and  t^  is  the  time  of  arrival  of  the  signal  at  microphone  i. 
If  we  select  tj  =  0  as  a  reference  and  use  six  microphones,  e  can  be  written 
as  follows. 

e  =  Y12(T12-(-t2))2  +  Y13(iir(-t3))2  +  •••  +  Y16(T16-(-t6))2 

6  6 

+  y  y  y. . [t. .-(t.-t.) ]2 

i>j  j=2  J  J  J 

A  necessary,  in  this  case  sufficient  because  of  convexity,  condition  for  e 
to  be  minimized  is  that 


i  =  2,3,. ..,6 


Taking  the  partial  derivitives  of  e  we  have  the  following  set  of  simultaneous 
linear  equations  in  t2,  t3,...,t  that  can  be  easily  solved. 


where 
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In  Che  present  edition  of  the  program  this  set  of  equations  was  solved  by 
using  a  canned  subroutine.  If  computer  storage  becomes  a  problem  the  solution 
vector'  can  be  obtained  by  using  a  form  of  the  steepest  descent  algorithm. 
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GLOSSARY  OF  PROGRAM  VARIABLES 


NMIC  =  Number  of  Microphones 

NL  =  Number  of  data  points  for  each  microphone 
NSW  =  Number  of  points  in  the  sliding  window 
NPFW  =  Number  of  points  in  the  Fourier  window 

SMIC(J,I)  =  Signal  amplitude  for  MIC  J  at  time  I  relative  to  1T(J) 

IT(I)  =  Starting  time  for  microphone  1  in  MSEC  relative  to  blast 
LRST(I)  =  Rough  starting  point  for  microphone  I  relative  to  IT(I) 

NR(I)  =  Starting  time  (MSEC)  for  Fourier  window  for  microphone  signal  I  rela¬ 
tive  to  IT(I) 

E(I)  =  Energy  within  correlation  window  for  signal  I 
TAU(J,I)  =  Raw  time  difference  between  windowed  signals  J  and  I 
RTAU(J,I)  =  Rough  time  difference  between  signal  from  microphones  J  and  I 
determined  by  correlation 

CTAU(J,I)  =  Corrected  time  difference  between  microphone  J  and  I  by  least 
squared  error  fit 

GAM(J,I)  =  Normalized  correlation  coef  between  signals  from  microphones  J  and  I 

G(J,I)  =  Reliability  weight  assigned  to  RTAU(J,I)  for  the  least  square  procedure 

AT(1)  =  Time  difference  shift  of  signal  I  relative  to  MIC  1  signal 

TIME(I)  =  Relative  time  of  the  Ith  signal  for  position  estimator 

A(J,I)  =  Augmented  coefficient  matrix  for  least  square  fit 

JC(I)  =  Control  variable  for  Subroutine  GJR 

V(I)  =  Variables  V(l)  and  V(2)  for  input  to  linear  equation  solution  subroutine  GJR 
CMEAS(I)  =  Overall  measure  of  correlation  for  microphone  I 


X(I)  =  Ith  value  of  first  signal  for  cross  correlation  on  input 
Ith  value  of  Fourier  transfer  on  output 
Y(I)  =  1th  value  of  second  signal  for  cross  correlation  on  input 
ith  value  of  Fourier  tranform  on  output 
WORK(I)  =  Work  space  specified  by  FOURC  subroutine 
XX(I)  =  Ith  value  of  estimate  of  noise  spectrum 
YY(I)  =  Ith  value  of  estimate  of  signal  plus  noise  spectrum 
ZZ(I)  =  Value  at  Ith  freq  for  prefilter 
IDAY(I)  =  Day  of  the  start  of  the  Ith  microphone  signal 
IHOUR(I)  =  Hour  of  the  start  of  the  Ith  microphone  signal 
IMIN(l)  =  Min  of  the  start  of  the  Ith  microphone  signal 
SEC(I)  =  Sec  of  the  start  of  the  Ith  microphone  signal 
MIC(I)  =  PASS  microphone  number  associated  with  mic  I 
UMAX  =  Microphone  number  selected  for  reference 
TEMP  =  Effective  temperature  at  time  of  shot  in  Deg.  C. 

MILS  =  The  effective  wind  direction  at  time  of  shot  in  mils 
KNOTS  =  The  effective  wind  speed  at  time  of  shot  in  knots 
KXX  =  Noise  window  starting  point  for  SMIC(I,L) 

KYY  =  Noise  window  starting  point  for  SMIC(J,L) 

PREF  =  Real  and  imaginary  value  of  prefilter 
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